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I. INTRODUCTION 



The free electron laser (PEL) was first proposed by J.M.J. Madey in 1971 [1] and lased 
in 1976 at Stanford University. The PEL technology was so promising that extensive funding 
and resources were devoted to PEL research as part of the Strategic Defense Initiative (SDI). 
Consequently, rapid developments in PEL theory and design were made throughout the 1980s. 
Currently, the technology has mamred to where high-average power PELs are feasible. 
Stanford University, Boeing, Rocketdyne, and the Continuous Electron Beam Accelerator 
Pacility (CEBAP) are among those pursuing kilowatt level PELs. 

The simplified schematic of a generic PEL is depicted in Pigure 1. 



bending 




mirror mirror 



Pigure 1 - Simplified PEL diagram 

A relativistic electron beam is deflected periodically by an oscillating magnetic field produced 
by a series of magnets called the undulator. The electron beam consists of a series of 
picosecond electron pulses, called micropulses, separated in time with a duty cycle about 10“^. 
Often the micropulses are produced in trains several microseconds long at repetition rates of a 
few hertz. One fuU cycle of this repetition is called a macropulse. However, continuous wave 
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(CW) operation of electron beams, where a continuous stream of micropuJses are injected into 
the undulator, can also be used. Light is radiated from the accelerated electrons, collected in an 
optical resonator cavity, and then amplified by subsequent electron pulses. 

Light recirculates within the optical resonator in pulses. The resonator cavity dimensions 
are aligned so that each optical pulse co-propagates with an electron micropulse as it travels 
down the length of the undulator. The micropulse is several hundred optical wavelengths long 
so that within a length of the micropulse corresponding to an optical wavelength, the electron 
longitudinal spatial distribution is essentially homogeneous. 

FELs are an attractive light source for several reasons. Perhaps the most positive 
attribute of FELs is. unlike chemical or solid-state lasers, that the optical wavelength can be 
tuned over a wide range of wavelengths. The wavelength is a funaion of electron beam 
energy. With slight variations in the beam energy, the FEL wavelength can be tuned over a 
range of about an order of magnitude about the center wavelength. A more complete 
description of FEL theory is presented in Chapter 11. 

The FEL also has disadvantages. The average power of FELs has been limited to only a 
few watts while the capital cost of the FEL, including the undulator magnets and a 
radiofrequency (rf) accelerator system used to produce and transport the electron beam, can 
cost tens of millions of dollars. Electron accelerators have been unable to produce the large 
average currents required to generate high-average power. Technology will now support a 
kilowatt scale FEL and a 100 kW FEL is expected within 5-10 years [2]. 

Free electron lasers are of interest to the Department of Defense due to their potential as 
high-power weapons and lower power countermeasure devices. A shipboard high-power FEL 
could provide essentially a zero time-of-flight defense against enemy missiles and aircraft. The 
need for such a weapon has been previously established [3][4]. Technological developments 
required to build an FEL weapon, however, are probably more than ten years away. 

A shipboard FEL would require an electrical prime power source. The amount and type 
of power required by an TOL would place a specific demand on a ship’s electrical power 
distribution system. Several possible choices of FEL system architectures are potential 
candidates for a high-power FEL. The maritime combat environment, however, places unique 
design constraints on the system. In Chapter III, the power requirements of the FEL 
architecture most suited to this environment are described. The methods used to determine 
these requirements are also valid for other architectures with appropriate changes in design 
parameters. All of the technologies considered are in use on existing FELs or have otherwise 
been demonstrated. No attempt to predict future technologies has been made. Chapter IV 
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outlines the prime power considerations in producing and distributing this power for the 
shipboard engineering plants that will power naval combatants into the 21®‘ century. 

The power analysis of Chapter III makes several assumptions that would affect the design 
and operation of the undulator and optical resonator. In Chapter V, simulations are presented 
to demonstrate that an FEL can operate at design levels under these assumptions. Current 
accelerator and optics technology, however, would not likely be able to support all of these 
parameters but it is important to validate the feasibility of pursuing such a design or a 
derivative. 

This thesis incorporates many of the aspects of FEL design. An FEL is not comprised of 
a set of unique characteristics. There are many design choices which must be made in order to 
define a given FEL architecture. The feasibility of a high-power shipboard FEL weapon is 
examined from a power analysis approach. A model for an FEL weapon is defined based on 
the volume, power distribution, and environmental constraints of the maritime combat 
environment. The power requirements for this model are determined and the characteristics of 
the required power distribution system are outlined. Finally, the optical characteristics and 
performance of this FEL model is explored. 
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II. FREE ELECTRON LASER THEORY 



The basic components of an PEL are a relativistic electron beam, a static periodic 
magnetic field, and a co-propagating optical field. The PEL operates through the quantum- 
mechanical processes of spontaneous emission, stimulated emission and stimulated absorption. 
The quantum-mechanical descriptions of these processes for bound electrons, as in an atomic 
laser, are well-known and easily formulated. It is not necessary, however, to use the 
quantum-mechanical descriptions for the PEL. The PEL can be completely described by 
classical electromagnetic theory relations since there is a high density of electrons within an 
PEL electron beam and photons within the optical wave [5]. The following theory has been 
collected from a variety of sources. The principal sources for each of the following sections 
are noted in the subject headings. The following relationships are presented in the cgs unit 
system. 



A. RESONANCE CONDITION [5][6] 

An accelerated electron radiates energy. In the PEL, the free electrons are accelerated by 
the static periodic magnetic field within the undulator. The electrons are traveling down the 
axis of the imdulator at relativistic speeds. Therefore, the period (wavelength) of the magnetic 
fields as they appear to the elecffon are not the same as the apparent undulator period as seen 
by an observer in the laboratory frame. This is often referred to as the relativistic Doppler 
shift. Since the undulator wavelength as seen by the electrons is a function of their energy, the 
wavelength of the photons emitted by the electrons is also a function of their energy. We can 
describe these effects, however, within the laboratory frame of reference due to the relativistic 
invariance of our classical relationships [7]. 

The electron beam interactions with the electromagnetic fields within the laser cavity are 
described by the Lorentz equations and the wave equation. The relativistic form of the Lorentz 
equations are; 

= -c(f -i-Kxfi) (2.1) 

at 

= -e t ( 2 . 2 ) 

di 



Y = 



(2.3) 



VTTf-f 

where y is the Lorentz factor, ^ = v^/c, v* is the electron velocity, c is the speed of light, m is 
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the rest mass of an electron, e is elementary charge magnitude of an electron, ^ is the electric 
field within the laser cavity, and ^ is the total of magnetic fields within the laser cavity. 

Equation 2.1 represents conservation of momentum for the system and is commonly 
known as the Lorentz force equation. Equation 2.2 desaibes energy transfer and would be 
called stimulated emission or absorption in a quantum-mechanical formulation. For an FEL we 
desire a net transfer of energy fi’om the electron beam to the optical field that is initially 
produced from spontaneous emission. This energy transfer is many orders of magnitude larger 
than the spontaneous emission contribution per pass of the electron beam through the 
undulator. We must ensure that stimulated emission occurs preferentially to stimulated 
absorption in order to make the FEL a useful device. 

The electric field, within the undulator is due to the presence of the optical field. 
Equation 2.2 shows that for ^ > 0, energy will be transferred from an electron to the optical 
field. Furthermore, special relativity states that the speed of an electron must be less than the 
speed of light in free space so that the electron must travel slower than ("slips" along) the co- 
propagating optical field. The transverse component of the electron velocity is periodic due to 
the acceleration of the electrons by the undulator. While traveling through one undulator 
wavelength, A.^, the electron’s transverse velocity, v^, will osciUate one full period as shown in 
Figure 2. In this figure, an electron’s longitudinal travel is shown with it’s average velocity 
over one undulator wavelength, PyC . After one undulator period, the electron slips behind the 
optical wave by one optical wavelength, X. At position z = 0, this electron is located where 
^ = 0 and V| = 0. As the electron moves down the undulator, increases and ^ becomes 
positive so that > 0 from 0 < 2 <Xol2. As the electron continues down the undulator, 
both and ^ are less than zero so that ^ ^ > 0 from Xol2 < z <Xo. The condition where 

the electron slips behind the optical wave by one optical wavelength as the electron traverses 
an undulator wavelength is called the resonance condition . A single resonant electron transfers 
energy from the beam to the optical field. 

The resonance condition defines an important feature of the FEL; specifically, its 
wavelength tunability. To illustrate this, we begin with the Lorentz factor equation. Since 
^ ^ = p.“ + the Lorentz faaor (equation 2.3) can be rearranged to solve for 

Pz = Vi "(1 +rPr)/r • (2.4) 

Define K~ = -y^p^ s undulator parameter. Typical values of K and y are ^ ~ 1 and y ~ 50- 
500 so that 1 -»■ y^PjV « 1. With a Taylor expansion, we have 
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X 
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Figure 2 - Electron slippage of an optical wavelengUi as it traverses a distance of 
one undulator wavelength. From Ref. [4] 



P, = 1 - . (2.5) 

2r 

The difference in velocity between the optical wave and the electrons is c(l-P, ). while the 
time it takes for the electron to traverse an undulator period is Therefore, the 

resonance condition is 






1 

2r 



(2.6) 



This mathematical representation of the resonance condition demonstrates the wavelength 
tunability of the FEL. In section II.C, we will see that 



K = 



2nmc' 



(2.7) 



where is the undulator rms magnetic field su^ength. The undulator field strengtli and 
wavelength characteristics are fixed after undulator consu^uction. Therefore, since X « 1/y^, tlie 
FEL optical wavelength can be rapidly adjusted over a wide range with small adjustments to 
the elecffon energy. Actual FEL experiments have produced optical radiation at wavelengths 
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of 10mm > >. > 300nm [5], 



B. SPONTANEOUS EMISSION [5][7] 

An accelerated electron will radiate electromagnetic energy. The power radiated by a 
relativistic electron is given by the generalized Larmor formula or Li6nard result. For an 
accelerating force perpendicular to the direction of electron motion, the total power radiated by 
a single electron is 

/’e (2.8) 



This power is distributed in space according to 

e _ 2e~Y^ j _ 4y^6~cos^(|> 

dQ. nc (1+7^02)3 (j^.y 202)2 



(2.9) 



where 6 is the angular deviation from the longitudinal axis in the plane of acceleration, ^ is the 
angular deviation from the longitudinal axis perpendicular to the plane of acceleration, and Q. is 
the solid angle into which the power is radiated. This power distribution is illustrated in Figure 
3. The total power is radiated into a narrow cone of width, A0 ~ 1/y and solid angle 
AD ~ 1/y^ centered along the direction of motion. This narrow beam causes the optical 
radiation within the FEL laser cavity to have an intense energy density complicating FEL 
resonator designs due to the limitations of the optical components. 

Photons emitted from the electron will have an energy, E = he IX, where X is given by 
the resonance condition. Therefore, 



E = 



2hck,y^ 

d+K^} 



( 2 . 10 ) 



where is the wavenumber associated with the undulator static magnetic field. The time the 
electron spends within the undulator interaction region is Ar = L/c, where L is the total length 
of the undulator. The number of photons emitted per electron in one pass of the undulator, 
then, can be easily determined by multiplying the total power emitted from an electron, 
equation 2.8, by At and then dividing by the energy per photon, equation 2.10. In the next 
section we will see that = K~koC~l-f so that the number of photons spontaneously emitted 
per pass, Wg , is 

Wg = ^ N oK^ (1+K-) 2N oK- (2.11) 
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normalized pattern 
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(a) (b) 

Figure 3 - Radiation pattern from a relativistic electron (y=100) under a 
transverse acceleration, (a) (j) = 0. (b) 4> = n/2 

where a = c^lTic » 1/137 is the fine suucture constant and N = L/X„ is the total number of 
undulator periods. For an undulator design with K ~ 1 and N = 100, about three photons will 
be emitted per elecu^on per pass through the undulator. Not all of these photons, however are 
emitted into the coherent optical mode. The coherent optical mode size is controlled by the 
Rayleigh length. z„, which is the propagation distance required for the optical mode cross- 
sectional area to double. The number of photons emitted into the solid angle AQ ~ MNy 
defined by the coherent optical mode with a Rayleigh length comparable to the length of tlie 
undulator [5] is 

VF,. = y o/rMl+A'^) « (2.12) 

Therefore, only about 1% of the electrons will spontaneously emit a photon into the coherent 
optical mode per pass. 

The above relations imply two significant factors. First, the total photon energy 
emissions are many orders of magnitude less than the total electron energy so that momentum 
recoil imparted to the elecU'on has a negligible effect on the elecu^on’s trajectory. Secondly, the 
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process of spontaneous emission constitutes the initial mechanism for FEL startup only. The 
primary power transfer mechanism from the electrons to the optical fields after a small power 
density has built up within the lasa resonator cavity is stimulated emission. At saturation, up 
to a few percent of the electron beam energy may be extracted per pass of a micropulse 
through the undulator by stimulated emission. 

The electron micropulses passing through the undulator create optical pulses of similar 
length. The expression for the resonance condition, equation 2.6, is derived based on the 
electron slipping one optical wavelength behind the photon every undulator wavelength. 
Therefore, we might expea that the minimum length of the optical pulse will be NX. This 
length is called the "sUppage distance" since the electron micropulse slips a total distance of 
NX along the optical pulse as the micropulse travels down the undulator. Acmal optical pulses 
are often longer than the shppage distance, and under certain conditions, called "short pulse 
effects," may even be shorter than the slippage distance. 

C. ELECTRON DYNAMICS AND THE PENDULUM EQUATION [5] 

The energy transfer between the electron beam and the optical wave occurs in the 
interaction region of the undulator. Within this region, the electric field is due exclusively to 
the optical electromagnetic wave while the magnetic fields are due both to the optical wave and 
the static undulator magnetic field. The first FEL was developed with helically polarized 
undulator magnetic fields and has a simpler mathematical devdopment for some FEL 
phenomena. Currently most FELs have linearly-polarized undulators. The fields within an 
undulator linearly polarized in the vertical (J) direction are illustrated by Figure 4 and are 
given by 



where (o = kc is the angular frequency of the optical radiation, the subsaipt "m " represents an 
undulator parameter and the subscript "r " represents an optical wave parameter. 

In order to describe the electron dynamics and trajectories, we return to the Lorentz 
equations as applied to a linearly polarized undulator. 



=^0 sm(k„z) y 
= E cos {kz -u>t -Ht)) X 
(kz -<at -h|)) y 



(2.13) 



(2.14) 



(2.15) 



at me L 



(2.16) 
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Figure 4 - Fundamental physics of the FEL interaction. After Ref. [4] 
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(2.17) 



vTT|=f 

Equations 2.16 and 2.18 combine to obtain 

Pi = cosik^z) 

y 

^ eBoK 

^ ~ 7 

2mnc 

Pz =Vl -tl +^^(1 +cos(2^„r)]/r 

and substituting equation 2.19 into equation 2.17 we obtain 

= V2 cos (kz -( 0 / +4) , 

dt yinc 



(2.18) 

(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 

( 2 . 22 ) 



Integrating equation 2.19 and plugging in representative values for existing FELs [5], we find 
that the transverse deflection of the electron beam as it passes through the undulator is on the 
order of 0.2 mm. Furthermore, imdulator wavelength values are on the order of 5 cm. 
Therefore, even with N = 100 periods for the undulator, the undulator is only a few meters 
long. The laser resonator, however, may be much longer to allow the optical beam size to 
increase and reduce the power density impacting on the mirrors. 



The electron trajectories described by the above equations have oscillatory components 
within an undulator wavelength. Over an undulator period, however, the electrons have an 
apparent average or constant trajectory. It is convenient to define the electron’s position in 
terms of the electron phase, as 



C = (/: + k^)z - (Of 



(2.23) 



where z is the position of the electron given its average motion neglecting the oscillations. 
Since k » ~ kAz. A change M ~ X corresponds to a A^ ~ 2it radians. From equation 

2.21, we obtain 



and 




(2.24) 




sin(.2ko ct) 



(2.25) 
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where ^ = A'^/2(l+A'*) and k^z ~ k^ct since the electrons are traveling at nearly the speed of 
light. Substituting equations 2.23 and 2.25 into equation 2.22 we obtain 



dt V2 yrnc 



cos[(C+4>)-^sin(2A:o ct )]+cos[(C+()))-2*o -^sin( 2 /:o ct )] 1. 






(2.26) 



Averaging equation 2.26 over one undulator period results in 



dt V 2 ynic 



cos (C-Ht>) 



(2.27) 



Differentiating equation 2.23 with respea to time twice and equation 2.24 once, they combine 
with the resonance condition and simplify to give 



d^-iQ _ ^k,eEK,(^) 
at Y'm 



COS(C-Ht)) 



(2.28) 



where K[^ (^) = K [J^ (^)-J i(^)]. Equation 2.28 is called the PEL pendulum equation . 

The pendulum equation describes an imponant physical phenomena in FELs called 
"electron bunching". This is easily seen in equation 2.27. Consider a group of electrons 
uniformly distributed over a = 2n space. Those electrons in the vicinity of ^+4> ~ n will 
experience a loss in energy and slow down. Those electrons near ^-hI) ~ 0 will gain energy 
and speed up. Therefore, the electrons will tend to bunch together near ^ = nil. Notice that 
the electron bunching is periodic over a = 2n corresponding to A: = The coupling 
between the pendulum equation and the wave equation allows for the net transfer of energy 
from a bunched electron beam. 



D. THE FEL SELF-CONSISTENT WAVE EQUATION [5][8] 



The propagation of the optical wave is governed by the wave equation. The form of the 
wave equation applicable in the undulator interaction area is 



C* dt- 



A = 



4n 7» 

“ ^ I 



(2.29) 



where A is the vector potential and T y is the transverse electron beam current. For an electron 
beam and optical wave in the interaaion region of a linearly polarized undulator 



^ ~ sin (k: -<ar +^) .x 


(2.30) 


Tyi = -ec'$y 6^(.v-r;) .v 


(2.31) 
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(2.32) 




cosik^z) 



X 



where 6 ^(a:) is the three-dimensional delta function and the subscript represents values of 
the electron such that / = 0 except at the exact location of an electron. The transverse 
beam current is a summation of the contributions of individual electrons. 



Substituting equation 2.30 into equation 2.29 results in a complex expression. Since we 
are using the PEL to produce coherent, monochromatic light, we can apply the slowly-varying 
amplitude and phase approximation. 



dE 


c k E , 


^ c- 

(o E , 
dr 




dz 




0(j) 

dz 


« A (f) , 


-e- 

3 

V 


(2.33) 



This approximation is useful for coherent monochromatic sources because the inverse Fourier 
transform of a narrowband signal corresponds to a slowly-varying envelope on the optical 
wave. This approximation simplifies the second-order wave equation to the first-order equation 



M+ 1 .^ 

dz c dt 



oos{kz -iol +^)-E 



1 

dz c dr 



sin(l:r -<D/+<t>) 




X . 



(2.34) 



Substituting in the expression for the electron beam current and averaging over a volume of the 
electron beam that is many optical wavelengths long, 2.34 simplifies to 



A lA 

dz c dr 



= -2<2enK 



^ -i la -i \(k-k, )z -^1 1 

Y, 






(2.35) 



where [ ~ ] represents the average over many optical wavelengths. Using equations 2.23 and 
2.25 we note that (l:-hA:o)r-(0/=^-^sin(21:o2) and (l:-A„)r-(i)r=^-2/:or-^sin(21:o- )• 
Computing the average in equation 2.35, the wave equation is now 



dz c dr 






Ee“^ = — b^{x-fj) . 

, Y- 



(2.36) 



An electron beam consists of a large number of electrons making the calculation of the 
summation in 2.36 prohibitive. However, we can take a "few" electrons over an optical 
wavelength representative of the entire electron beam within that wavelength and compute the 
summation. Therefore, the FEL wave equation is 



A+li. 

dz c dr 



Ee'* = -l^neKLi^) 



< p(=)<’"'^/Y>,=z-p^„ 



(2.37) 
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1 ^ 

where p is the electron density and <...> = — 52(...). 

P /=i 

The wave equation and the pendulum equation are the two important equations describing 
the electron beam and optical wave interactions. However, it is more convenient to non- 
dimensionalize these equations for PEL analysis. 

E. NON-DIMENSIONAL ANALYSIS AND PHASE-SPACE [5] 

The pendulum equation, 2.28, and the wave equation, 2.37, are the principal equations of 
motion describing the electron beam and optical wave dynamics within an PEL. In most 
instances, it is necessary to analyze PEL properties numerically. Numerical analysis is simpler 
and results more physically meaningful when these equations are used in a non-dimensional 
form. 

A convenient choice for non-dimensionalizing the PEL equations is to establish a 
coordinate system by normalizing time to the time it takes the optical wave (or electron beam) 
to traverse the length of the undulator and to reference distance according to the position of a 
point traveling at the speed of light. The transformations are t = c/L i and : = ~ - c t . where 
X is the non-dimensional time and 5 is the position with respect to the reference position. The 
coordinate 2 stays with a point on the field envelope traveling at speed c. Por the special case 
where the electron pulse is many optical wavelengths long so that p(r ) = p and the electrons 
have a uniform and constant energy, then the non-dimensional pendulum equation and wave 
equation become 



a’-r l^lnU^eEK^i^) 

TTt = C = , , CO.V (C-Hj)) 

dT“ 


(2.38) 


— ^{Ee'^) = 2<2neKL{^)^ <£>-'?> 
L ax Y 


(2.39) 



Introducing the complex optical field envelope a = tr; I and the dimensionless 
current, j, where \a\ = 2^nLNeEKi{^)l-fmc~ and j = the non- 

dimensional pendulum equation and wave equation are 

OO 

^ = la I cos(C-Ht») (2.40) 

a = - j < c"? > . (2.41) 

Typically, the optical mode waist (radius). w„ = slkz^ln, within the undulator is larger than the 
radius of the electron beam. Therefore, only pan of the optical wave is amplified by the 
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electrons. This is compensated for by including the "filling factor," F = in the 

dimensionless current, j ^ jF . 

The optical wave and electron evolution within the electron beam are coupled according 
to equations 2.40 and 2.41. The electron phase, C. represents the electron position within the 
electron pulse and the change in electron phase, ~ A:Ar , is indicative of the change in 
electron position. The electron phase is periodic over a distance corresponding to one optical 
wavelength. Therefore, if we analyze the C evolution over a length of one optical wavelength 
(2n radians in phase) of the electron beam, then we know the C evolution over the entire 
electron pulse. It is convenient to select the ^ reference so that <() = 0. 

Another important parameter in our analysis of electron evolution is the phase velocity, 

O _ 

V = ^ = L[(k kg)^. - ^]. Differentiating equation 2.24 and non-dimensionalizing we find 
that Av = 4ji/V Ay/Y- Therefore, the phase velocity is indicative of the electron energy and the 
change in phase velocity is indicative of the change in electron energy. A phase velocity of 

V = 0 corresponds to an electron exactly on resonance, whereas an electron with a v = rt will 
advance about XU ahead of an electron on resonance over the length of the undulator. 

The FEL phase space evolution depicted in Figure 5 shows an example of the change in 
electron energy and position for an initial uniformly distributed monoenergetic electron beam. 
If we can form electron bunching around ^ = n, then a net amount of energy wifi be 
transferred from the electron beam to the optical field. An electron bunch centered around 
C = 0, however, will cause a net energy to be absorbed by the beam from the field. Figure 5 
shows the phase space evolution of a resonant electron beam, and a beam slightly above 
resonance. A resonant electron beam has no gain. This is because for a uniformly distributed 
electron beam at resonance, every electron which contributes energy to the optical field is 
balanced by another electron which absorbs an equal amount of energy from the field. An 
electron slightly above resonance will contribute a net energy to the optical field. This will be 
explained more fully in the next section. Through the process of mode competition in the 
FEL, the optical wavelength will adjust itself so that the electron beam will be slightly above 
resonance. 

The "football-shaped" curve in Figure 5 is called the separatrix. The separatrix curve is 
given by 

v/ = 2\a\ [1 - since, + ^)] (2.42) 

where the subscript "s " refers to the separatrix values. An electron orbit enclosed within the 
separatrix is a "closed" orbit. This is equivalent to a pendulum swinging side-to-side with 
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Figure 5 - Phase space evolution of an electron pulse, (a) resonance, (b) maximum 
gain 
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insufficient energy to reach the vertical position and "swing over the top". An electron orbit 
outside of the separatrix is an "open" orbit. Since these changes are periodic, an electron 
which leaves the field of = 2n shown is replaced by another electron from an adjacent 
A^ = 2n section of the micropulse. This is equivalent to a pendulum swinging over the top on 
every revolution due to a total energy greater than that required to overcome the potential 
energy of the vertical position. In our FEL, however, energy is being transferred both to and 
from the electrons so they are not necessarily confined exclusively to one region or the other. 
As the complex field strength, \a I, increases, the separatrix becomes larger trapping a larger 
fraaion of the electrons. For a maximum tfansfer of energy from the electron beam to the 
optical field, the electrons must bunch at ^ = n in the lower right-hand section of phase space 
enclosed within the separatrix, 

F. WEAK OPTICAL FIELDS EVOLUTION [5] 

From the pendulum and wave equations, it is clear that electron bunching around ^-Ht) » n 
provides for maximum gain of the optical fields. However, if the electron beam is over- 
bunched such that the bunch moves toward ^+<t) ss 0 in phase space, then energy will be 
absorbed back into the electron beam. The weak-fields assumption, la I $ n, is defined such 
that the electron beam is not over-bunched. Over-bunching is indicative of strong optical field 
FEL interactions. 

Under weak fields, we can expand the pendulum equation and wave equation to first 
order in la I and combine to form 

a = -^ f ri F(t') la(T-ri)l e dx' (2.43) 

2 Jo 

where 

oo 

F(x')= dq. (2.44) 

Equation 2.43 is referred to as the integral equation . F{x') is the characteristic function of the 
beam quality distribution f (q). The initial phase velocities of the electrons are v, = + q 

and are distributed about according to / (q ). For a perfect beam, where all of the electrons 
are injected monoenergetically with phase velocity into the undulator directly on axis with 
no transverse position or velocity components, then F (x) = 1 . 

The wave equation, 2.41, tells us that for low current conditions, j £ n, the optical field 
does not change significantly per pass. Assuming perfect beam quality, the solution to the 
integral equation gives an optical field development per pass to first order in j , 
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(2.45) 



a ( t ) = ci^\ \ + j 



2 - 2cos(Vq t) - v„ Tsin(v„ x) 



and 



(!)(x) = j 



2sin(VoX) - VqT(1 + cos(VqX)) 

2v.^ 



(2.46) 



where is the initial complex field strength. The gain of the optical field through the 
undulator in this low gain region is given by 



G (v„ ,x) 




9 





2-2cos (v„ x)-Vo isin (v^ x) . 



(2.47) 



The weak fields, low current gain spectrum is shown in Figure 6. 



0.135 

0.0 

- 0.135 
0.081 

0.0 

- 0.081 

Figure 6 - Weak fields, low current anti-symmetric gain spectnim 




The anti-symmetric gain curve shows several important FEL properties. First, there is no gain 
at resonance, Vq = 0. Secondly, the electron phase velocity must be slightly above resonance 
to achieve gain. The maximum gain is C = 0.135v and occurs at an initial phase velocity v^, = 
2 . 6 . 
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Under high current, high gain conditions, j » n, it can be seen from the wave equation, 
2.41, that the magnitude of the optical field changes appreciably over a single pass. The self- 
consistent solution to the integral equation for this region provides much different results. 
Maintaining the assumption of perfect beam quality, an PEL operating on resonance, = 0, 
will have a complex field development given by 



3 



<Kt) ~ 






± 

2 



2 



giving gain along the undulator of 

G(t) = 1 ^072)''^'3x 



(2.48) 

(2.49) 



(2.50) 



There is now gain for an electron beam directly on resonance. Off resonance, the single-pass 
gain at the end of the undulator is given by 



G 




072)''V3(1 -(v„t«/3)=) 



(2.51) 



where Tg = (2/y)'^^ is the bunching time. The gain spectrum for the high gain region is shown 
in Figure 7. As in the low gain case, the peak gain is still slightly above resonance. The phase 
velocity that gives the peak gain, however, is now a function of current. As j — > the gain 
spectrum approaches a symmetric spectrum with the peak gain occuring at resonance. 

The gain spectrum bandwidth can be determined from equation 2.51. Defining the 
bandwidth as the range of phase velocities where the growth rate is reduced by 10% [5], the 
bandwidth is Av„ = Using this same definition, the gain bandwidth for the low gain case 
is about n/2. Therefore, the bandwidth for the high gain case is much larger than for the low' 
gain case. 



G. GAIN DEGRADATION DUE TO BEAM QUALITY [5] 

The previous gain analysis assumed perfect beam quality. A real electron beam has an 
emittance, e, where the electrons enter the undulator at an angle. 0j, with respect to the 
undulator axis, or at a transverse position slightly off axis, . The elecnons will also have a 
small energy spread between them. The emittance and energy spread of the electron beam 
translates into a spread in phase velocities in phase-space. The phase velocity spread of an 
imperfect electron beam due to emittance is given by 
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Figure 7 - Weak fields, high current gain spectrum 



. 27TTV O'). 

~ "T^ V + r0.v') 
i+/C“ 



(2.52) 



and the phase velocity spread due to energy spread is Av = 471^ Ay/y. 

Imperfect beam quality reduces the single-pass FEL gain by degrading electron bunching. 
Consider the integrand of the integral equation. The complex field strength grows as the 
electron beam travels down the undulator causing the integrand to increase. Therefore, the rate 
at which the optical field strength grows is faster. Equation 2.51 and a plot of equation 2.47 
shows that this feedback causes exponential growth in the FEL optical field strength. With 
imperfect beam quality, however, the F(x) decays over time so that this feedback mechanism is 
attenuated. 

Electrons enter the undulator with a spread in phase velocities due to a spread in 
energies. Typical energy spreads from an accelerator are on the order of Ay/y « 0.001. The 
typical distribution of electron phase velocities is Gaussian about , 

-rn-oh 

fc(q) = -^ — (2.53) 

y2n 0 (j 
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where Oq = 4jiN Ayly is the standard deviation. The corresponding characteristic function is 

FciX) = e"''^^'-. (2.54) 

Fq (t) decays over the length of the undulator with a time on the order of 1 /oc . 

The change in phase velocity due to beam emittance has a quadratic dependence on 
electron position and entrance angle. Therefore, f {q) is not the same as the distribution of 
and 0j. Electron pulses normally have a Gaussian distribution in space. A beam with matched 
contributions from position and angle. = KkgX„, where 6 ^ and are the rms deviations of 
the electron beam entrance angle and position, will have an exponential phase velocity 
distribution given by 

= <7<0 (2.55) 

Oo 

where /o(<?) = 0 for <7 >0 and Oq = The corresponding characteristic 

function is 

Fq{-^) = ^. (2.56) 

F q{x) decays over the length of the undulator with a time on the order of I/O 0 . A Oq or Oq of 
about n will cause a random phase spread of about n. 

The effect of beam quality on gain is illustrated in Figure 8 for the case of electron 
energy spread. At Oq = 0, the gain spectrum is identical to Figure 6 . As Oq increases, 
however, the gain spectrum broadens, the phase velocity for peak gain increases, and most 
importantly the peak gain is reduced. Electron beam emittance has a similar effects. 

H. STRONG OPTICAL FIELDS EVOLUTION AND SATURATION [5] 

When an FEL reaches the strong field regime, the gain per pass is reduced. The FEL 
reaches saturation when the gain decreases to where it equals the optical resonator loss rate 
including the outcoupling of the mirrors. Strong fields can be detected in phase space by over- 
bunching of the electrons. This occurs when Ir; > n. 

In strong fields, the separatrix now traps a significant amount of electrons in closed 
orbits. An electron bunch is quickly formed around ^ + (J) = n where these electrons amplify 
the optical wave. As the electrons lose energy, their phase slips towards the ^ -t <t> = 0 position 
and energy is absorbed back from the optical wave. This reduces the net gain per pass of the 
FEL. 
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Figarc 8 - Degradation of the weak fields, low current anti-symmetric gain 
spectrum due to electron beam energy spread 

The gain degradation due to strong fields is illustrated in Figure 9. As the complex field 
strength and height of the separatrix increases, so does the gain degradation. Careful 
examination of Figure 9 reveals that strong-fields produce similar gain degradation 
characteristics as the energy spread effects of Figure 8. Tliis is because the FEL interaction 
induces an energy spread among the electrons. Therefore, the gain is also reduced in a manner 
similar to a phase velocity spread of Av = 4nN Ay/y. As the total energy loss of the electrons 
increases, a larger phase velocity spread is induced on the electrons. Eventually many electrons 
will move near the ends or even outside of the gain spectrum bandwidtli reducing the gain of 
the optical fields. As the optical fields approach saturation, more total energy is still being 
removed from each passing micropulse, but the fractional increase in optical energy (gain) is 
being reduced. 



23 



1.1 




Figure 9 - Gain degradation due to strong optical fields 



24 



III. FEL SYSTEM POWER REQUIREMENTS 



The principles of FEL operation outlined in the previous chapter describe how energy is 
extracted from a relativistic electron beam and radiated as coherent optical light. The 
production of a high quality relativistic electron beam, however, is a necessary engineering 
requirement. The basic functional components for a high-power FEL are the electron injector, 
the linear accelerator, the undulator and optical resonator, and the electron beam dump. Many 
technologies and configurations for each of these elements have been developed and applied. 
Existing FELs, however, operate at only a few watts of optical output and do not suffer from 
the infrastructure constraints of the naval shipboard environment including space availability, 
power generation and distribution, and adverse environmental conditions. In this chapter, the 
power requirements for a high-power FEL as applied to a potential shipboard weapons system 
are explored. All of the following technologies are currently in use or under development for 
FEL applications. Some of these technologies, however, will need further engineering 
development to support high-power operation. 

A. SHIPBOARD ARCHITECTURE 

In developing a shipboard high-power FEL, the issues that are of principal concern are 
system power efficiency, system size and weight, and personnel radiation hazards. The system 
size and weight are a direct function of beam energy, linear accelerator rf field frequency, 
system power efficiency, and optical resonator power density. Electron beam energy and the rf 
frequency determine the linear accelerator size and weight while the system power efficiency 
and the rf frequency determine the rf power source size and weight characteristics. Klystrons, 
klystrodes, and other high-power rf power sources suitable for accelerator applications 
contribute a significant fraction of the system size and weight. The optical resonator power 
density is much higher in an FEL than for a chemical or solid-state laser of equal energy 
storage due to the very small beam size. Therefore, FEL resonator length tends to be large in 
order to reduce the power density at the mirrors. 

Of the existing or proposed, FEL architectures, only the energy recovery system depicted 
in Figure 10 optimizes all three concerns. In this scenario, electron pulses are injected into the 
accelerator at some small initial energy of approximately 4 MeV. The accelerator increases the 
beam energy by an order of magnitude or more. Bending magnets guide the beam and direct it 
into the undulator. The undulator wiD remove approximately 2 percent of the beam energy and 
convert it into optical energy. Additional bending magnets then guide the beam back into the 
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optical 




beam dump 



Figure 10 - Generic energy recovery FEL arcliitecture 



accelerator where it enters 180 degrees out of phase with the accelerating fields. The 
decelerating electron beam then transfers it’s energy back into the rf fields. The energy 
recovered from the decelerated electron beam is then used to accelerate another electron pulse 
to repeat the cycle. The decelerated electron pulses are then directed into a beam dump where 
their residual energy is dissipated. This residual energy is about 5 MeV and is below the 
threshold of generating neutron radiation when it is dissipated in the beam dump. 

This architecture is the most appropriate for several reasons. First, the energy efficiency 
is increased dramatically and the physical size of the system is reduced to that required in a 
shipboard environment. Most importantly, it is the only architecture which would reduce the 
personnel radiation hazards to manageable levels since there will be no neutrons and much 
lower energy gamma radiation. Table 1 is a listing of the neutron production threshold for 
several elements. These elements are commonly used as structural materials in shipboard 
applications or have significance with respect to neutron production in accelerator beam dumps. 
Only beryllium (Be) has a neutron production threshold below 5.9 MeV. Therefore, as long as 
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Element 


Threshold 

(MeV) 


Element 


Threshold 

(MeV) 


Element 


Threshold 

(MeV) 


Be 


1.7 


Cu 


9.9 


Mo 


7.4 


Os 


5.9 


Ni 


11.4 


Mn 


10.2 


*Li 


7.3 


Fe 


11.2 


Co 


10.5 


*C 


18.7 


Zn 


10.2 


Cr 


12.0 



Table 1 - neutron production thresholds [9]. All elements that have isotopes of 
> 10% abundance with thresholds below 6 MeV are included. * denotes 
elements with isotopes of less than 10% abundance with thresholds below 
6 MeV 



beryllium is excluded from the beam dump structural materials, neutron generation is not an 
issue. An architecture which does not utilize energy recovery will contain an inherent large 
radiation problem that will be much more difficult to solve due to the large radiation energies 
involved. 

The feedback stability of this type of system has been numerically modeled and shown to 
be stable for several designs [10]. However, more comprehensive simulations, experiments, 
and possibly control systems engineering will be required to ensure stability at the higher 
currents examined in this thesis. 

Since we are concerned with the FEL’s potential as a naval laser weapon, the system 
must produce high output power at a wavelength suitable for propagation in a maritime 
atmospheric environment. The following power analysis will be developed for an FEL with a 1 
MW optical output power at a wavelength of 3 pm. The FEL wavelength can be tuned over a 
range from this operating point so we will stUl have some flexibility in the actual operating 
wavelength. To establish the operational requirements, it is assumed that this laser will only 
need to operate for three shots of 15 second duration within one minute and ten shots over a 
twenty minute interval. This requirement has no basis in the time required to destroy a missile 
with a laser. Rather, it was chosen because a missile U'aveling at a speed of 0.9 Mach could be 
engaged over a closing distance of about 5 km. 

To estimate the system power requirements, it is necessary to determine the electron 
beam energy and current. The resonance condition of equation 2.6 gives the required 
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relationship between the electron beam energy, the undulator spatial wavelength, and the 
optical wavelength. For given undulator and optical wavelengths, the electron beam Lorentz 
factor must be 



Current and planned FELs are summarized in Table 2 [11], The undulator parameter, K, is 
typically on the order of unity so that y » 

The power in an electron beam is the product of the electron beam energy, = ytnc^, 
and the average electron beam current, lavg. For the architecture of Figure 10, only about 2 
percent of the beam energy can be extracted by the undulator to allow for efficient energy 
recovery [3] [12] [13]. However, only a few experiments have been performed which 
demonstrate energy recovery but none of these have determined the accelerator energy recovery 
efficiency, 11^^^ , dependence on the undulator energy extraction efficiency. Ti^, . Therefore, the 
assumption that no more than 2% of the beam energy can be extracted within an energy 
recovery architecture may be conservative. Regardless of t|o, an elecu-on beam cannot be 
decelerated much below 3 MeV because the electron velocity and rf phase velocity are no 
longer essentially the same. The maximum energy recovery efficiency is given by 
Tlma,x = 1 - (3 Mev )IW ^ . However, energy extraction in the wiggler introduces an energy 
(velocity) spread among the electrons in the beam which further reduces y\rec- More 
experiments are needed to determine the dependence of on Tio . 

For a 1 MW optical output with tip = 2%, 50 MW of electron beam power is required. 
In order to limit the undulator and optical resonator size and the required beam current, a 
maximum beam energy is desired. From equation 3.1, Wj, = 100 MeV and la^g = 500 mA are 
chosen as the design parameters for the production of 50 MW of beam power. For this beam 
energy, a = 95% is assumed. 

It is desirable to maximize so that the total electron beam power can be reduced. 
This will provide several advantages. First, by reducing the required lavg , a high-power FEL is 
closer to realization. Current technologies can only expect to produce average currents of 10- 
100 mA [2][14]. Second, reducing reduces the size, weight, and expense of the accelerator 
structures and allows for the production of longer wavelength light that may be more 
advantageous in a maritime environment. 




(3.1) 
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FEL 


X 

(Rm) 


O: 


Ea 

(MeV) 


/ 

(A) 


N 


(cm) 


K 


EXISTING FELs: 
















UCSB(mm FEL) 


338 


25ps 


6 


2 


42 


7.1 


0.71 


S tanford(FIRFEL) 


80-200 


15ps 


4 


8 


50 


1 


0.7 


UCSB(FIR FEL) 


63 


25ps 


6 


2 


150 


2.0 


0.13 


Tdcy o(TJT-FEL) 


43 


lOps 


13 


22 


40 


4 


0.7 


Holland(FELIX) 


40 


5ps 


25 


50 


38 


6.5 


1.5 


Osaka(ISIR) 


40 


30ps 


17 


50 


32 


6 


1 


Bruyeres(ELSA) 


20 


3(^s 


18 


100 


30 


3.2 


0.8 


Frascati (LISA) 


15 


7ps 


25 


5 


50 


4.4 


1 


Grumman(CIRITD 


14 


5ps 


14 


150 


73 


1.36 


0.2 


Beijing(IHEP) 


10 


4ps 


30 


14 


50 


3 


1 


Orsay(CXIO) 


8 


0.3ps 


50 


80 


48 


4 


1 


LANL(AFEL) 


4-6 


lOps 


15 


200 


24 


1 


0.27 


Darmstadt(IR-FEL) 


5 


2ps 


40 


2.7 


80 


3.2 


1 


Stanford(SCAFEL) 


5 


0.7ps 


37 


10 


72 


3.1 


0.83 


Vanderbilt(FELI) 


3 


Ips 


43 


50 


47 


2.3 


1 


EhikeCMarkEH) 


3 


3ps 


44 


20 


41 


2.3 


1 


BNL(ATF) 


0.5 


6ps 


50 


100 


70 


0.88 


0.35 


BNL(ATF-UV) 


0.25 


6ps 


70 


100 


70 


0.88 


0.35 


LANL(APEX) 


0.37 


lOps 


46 


135 


73 


1.36 


0.58 


Tsukuba(NUl-IV) 


0.35 


160ps 


300 


5 


2x42 


7.2 


2 


Orsay(Super-ACO) 


0.35 


250ps 


800 


.1 


2X10 


13 


4 


Okazaki (UVSOR) 


0.3 


6ps 


500 


5 


2x8 


11 


2 


Novosibirsk(VEPP) 


0.24 


35ps 


350 


6 


2x33 


10 


1.6 


PROPOSED FELs: 
















Florida(QREOL) 


231-600 


CW 


1.7 


0.2 


185 


0.8 


.15 


Netherlands(TEUFEL) 


180 


20ps 


6 


350 


50 


2.5 


1 


Rutgers(IRFEL) 


140 


25ps 


38 


1.4 


50 


20 


1 


Moscow (Lebedev) 


100 


20ps 


30 


0.25 


35 


3.2 


0.75 


Tokai(SCARLET) 


40 


4(^s 


15 


10 


62 


3.3 


1 


Stanford(FIREFLY) 


40 


2ps 


20 


6 


25 


6 


1 


LBL(IRFEL) 


3-50 


33ps 


55 


60 


40 


5 


1 


CEBAF(IRFEL) 


2-24 


2ps 


50 


100 


25 


6 


3.1 


Boeing(APLE) 


10 


60ps 


17 


140 


101 


2.4 


0.2 


Boeing(APLE) 


10 


18ps 


34 


450 


257 


3.9 


1.2 


Stanford(FEL) 


10 


4ps 


24 


25 


52 


2.6 


0.87 


Osaka(ILT) 


10 


2ps 


9 


100 


30 


.66 


0.3 


UCLA(IRFEL) 


10 


2ps 


20 


200 


40 


1.5 


1 


Novosibirsk(RTM) 


7 


50ps 


51 


100 


4x40 


9 


1-2 


BNL(HGHG) 


3.35 


lOps 


30 


110 


83 


1.8 


1.4 


CEBAFdJVFEL) 


0.15-2 


0.4ps 


200 


200 


48 


3 


1.5 


Osaka(FELI) 


1 


2ps 


170 


100 


50 


6 


1.26 


Rocketdyne 


0.84 


3ps 


90 


500 


160 


2.4 


1.4 


Dortmund(DELTA) 


0.4 


5(^s 


500 


90 


17 


25 


3.1 


ENEA-Frascati 


0.2-0.3 


15ps 


2.3 


4 


8 


2.5 


1 


Haiima(HID 


0.28 


lOOps 


500 


3 


170 


1.8 


4.2 


BNL(DUVFEL) 


0.075 


6ps 


310 


300 


682 


2.2 


1.54 


Duke(Ring) 


0.05 


lOps 


1000 


350 


2x33 


10 


1.7 


SLAC(L(XS) 


0.0004 


O.lps 


7000 


2500 


723 


8.3 


4.4 1 



Table 2 - Summary of existing and proposed FELs [11] 
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B. ACCELERATOR STRUCTURE 

A typical linear accelerator is comprised of a number of rf cavity cells like those depicted 
in Figure 11. 




B-Faclorv Cell 
Shape 




Copper 
Storage Rino 
C.gll Shagg 



(a) (b) 

Figure 11 - Typical accelerator cells, (a) superconducting (b) room temperature. 

From Ref. [15] 

Linear accelerator structures are fabricated from both room temperature (RT) and 
superconducting (SC) materials. SC structures reduce operating power requirements due to the 
low ohmic losses in the accelerating cavity walls. They also have higher allowed operating 
fields and accelerating gradients as well as much larger apertures between cells in the cavities. 
An accelerating electron forms an image charge within the accelerator structure. As it passes 
through an aperture connecting two cells, the image charge is now much closer to the electron 
beam axis. These image charges have a disruptive effect on the beam quality. The larger 
apertures mean that image charges will have a lesser disruptive effect in SC structures than RT 
structures. The higher gradients mean that SC structures can accelerate the electron beam to 
the desired W,, over shorter accelerator lengths than RT structures reducing the total weight 
and volume of the structure. RT structures, however, are forced to have geometries that 
intensify the harmful beam quality effects in an effort to minimize power dissipation [16]. RT 
structures still have applications because they are less expensive and easier to fabricate than the 
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sc structures and they don’t require a large and inefficient liquid helium refrigeration system. 



The development, theory, and design behind accelerators is complicated and requires 
extensive numerical modeling. In performing a power requirements analysis, however, 
accelerating cavities can be characterized by a few parameters for most purposes [16] [17] [18]. 
The significant parameters are the rf frequency, / (MHz), the cavity quality factor, Q , the 
shunt impedance, Z = IQ (Q.lm ) where R^h is the shunt resistance, the cavity active length. 
La (m), and the average acceleration gradient. £<, (MV/m). The power dissipated per unit 
length is 



P = 



Z Q 



(3.2) 



RT cavities typically have a 0 ~ 2x10’’ and an == 2 MV/m compared to values for SC 
cavities of 3x10*^ and 10 MV/m respectively. Some actual values are shown in Table 3. The 
dissipated power is supplied to the accelerator by the rf power source. 



cavity 


/ 

(MHz) 


Ea Q 

(MV/m) 


Z 

(Qlm) 


TESLA (SC) [19] 


1300 


20 


>5x10’ 


1100 


CEBAF (SC) [2] 


1500 


8 


5x10’ 


1020 


CDRL (SO* [20] 


500 


5.3 


2x10’ 


400 


TEUFEL (RT) [21] 


1300 


2 


18330 


1714 



Table 3 - Typical cavity parameters (* - parameters from a proposal) 



In order to properly apply equation 3.2, it is necessaiy to more fully characterize the 
parameters. Superconductivity is charaaerized by the BCS (Bardeen, Cooper, and Schrieffer) 
theory. For a SC cavity, the surface resistance, R^ = Rbcs given by the relation 

[16][17] 

A(n 

Rs=A + /?,,, (3.3) 

where A (Q K/m/MHz-) is the BCS material-dependent constant, A (eV) is the energy gap of 
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the SC material used to fonn Cooper pairs, T (K) is the absolute temperature of the SC 
structure, kg (eV/K) is the Boltzmann constant, and is the residual losses of the SC 

structure after fabrication ( ~ 1-10 nCl/m ). The quality factor of an accelerator cavity is 

the ratio of the energy stored within the cavity to the if energy dissipation within the cavity per 
rf cycle. The Q value is related to the surface resistance through a geometry factor, G , by 




Q QbCS Qres 



(3.4) 



where Qscs ‘s the contribution to the quality faaor due to losses in the superconductor and 
Qres is the residual losses contribution to the quality faaor. The geometry factor varies from 
270 to 290 ohms [17]. An average value of G = 280 £2 is assumed. With proper fabrication 
and conditioning techniques, a Qrgs, of better than 5xl0’ can be obtained [16][17]. 

The shunt impedance is also a function of frequency. In the absence of other changes, 
the shunt impedance is proportional to frequency [22]. However, as the frequency changes, the 
subsequent scaling of the cavity dimensions allow us to make more favorable, slight 
adjustments in cavity geometry. A curve fit to existing SC cavities is [17] 

Z =380(7/500)'’ (O/m). (3.5) 



The cavity geometry changes do not have a very large effea on the shunt impedance. This is 
due to charaaeristics of SC cavities like the very large apertures which can be seen in Figure 
11. As a result, the variations are incorporated into the frequency scaling of equation 3.5 
instead of the geometry factor. 

SC cavities are currently made exclusively of niobium (Nb). Although there are many 
SC materials available, it is desirable to use a material with the highest superconducting critical 
temperature, (Tg), and that can support the highest rf field strength, (Hg), without forcing it out 
of its superconducting state. It must also have a widespread availability, relative ease of 
fabrication, high thermal conductivity to dissipate heat, significant material phase stability, and 
have stable superconducting properties. Only lead (Pb) and Nb currently meet these 
requirements with Nb having the larger values of Tg and Hg. High temperature 
superconductors are currently unsuitable for these applications because of large surface 
resistances (Rs ^ lO'^Q) and low Hg (the highest values are 1/16 that of Nb) [16]. 

The BCS material constant for Nb is /4 = 00(X)28/(5()0)‘. The static heat load of a SC 
cavity due to power couplers, etc. is typically about 8 W/m at 500 MHz [17]. Therefore, the 
total power dissipated in the cavity is = P static + ^bcs + f^res < that is 
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£/(/■ /500)‘'(.000028)f>-‘’ 
380(280)7 



(W/m) . (3.6) 



RT cavities, normally fabricated from copper, have a surface resistance which is 
proportional to the square root of the frequency and the cavity Q is inversely proponional to 
the square root of the frequency for a fixed geometry [22] [23]. Therefore, 



where I'l and A: 2 ^re geometry constants that are inversely proportional to the square root of 
both the electrical resistivity and magnetic permeability of the cavity material. The shunt 
impedance, then, is proportional to frequency and independent of the material from which the 
cavity is fabricated or its surface finish [22]. The constants, however, are heavily dependent on 
cavity geometry. Unlike with SC structures, the geometry constants vary widely with the rf 
frequency as adjustments to cavity geometry are made with the scaling of the cavity 
dimensions. This is apparent from Figure 11 with the RT cavity cell’s narrow apertures. The 
values of for existing designs vary over the range 0.9 < A'l < 4 (12/m/MHz *^") 
[12] [20] [21] [23]. Lower frequency cavities will have larger values of Ati. 

The electron beam is a series of micropulses that are separated in space by some integer 
number of rf wavelengths, X^f . These micropulses can either be injeaed into the accelerator in 
macropulses or in CW mode. For a high-power FEL, CW mode is preferred [13] because: 

1. The peak power of a micropulse is reduced. In reducing the total charge contained 
within a micropulse, it is easier to overcome the disruptive space charge effects that 
degrade the electron beam quality. Furthermore, the problem of electron beam "halo" is 
reduced. A real electron beam has a uansverse spatial disuibution. Although the electron 
density away from the beam transpon axis is a small fraction of that on axis, the number 
of these electrons in a high current pulse is significant. These electrons can strike the 
magnetic beam transport optics and deflect so that they excite rf modes and cause 
transport instabilities for the entire beam. Similarly, the activation of the optics would 
create a significant radiation hazard over a very short time [12][14]. 




(3.7) 



and 




(3.8) 
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2. The peak power of the optical pulses developed in the optical resonator is reduced 
extending the lifetime of the mirrors and various optical components. Furthermore, the 
more constant loading of the optics reduces the need for stabilization measures due to 
transient deformations. 

3. The overall power efficiency of the accelerator is increased. When electron pulses are 
accelerated in macropulses, the rf power is applied to the accelerator in pulses (pulsed 
power operation). A finite "filling time" is required to build up the energy storage within 
the accelerator cavity prior to the injection of the next macropulse. There is significant 
energy lost in the accelerator during both this filling time and the decay time at the end 
of the macropulse where stored energy within the cavity is dissipated. 

4. Optimization of an energy recovery system is possible. Pulsed power systems have 
significant time variations of the electron beam parameters over the duration of the 
macropulse. This affects the optical evolution in the resonator. In CW mode, these 
variations occur over only a small firaction of the electron beam’s duration and are not 
repetitive. 

5. Finally, the control systems are easier to realize. They can be simpler, slower, and much 
less expensive than for an equivalent pulsed power machine. This is an extremely 
important factor because the feedback and control requirements for a high-average power 
energy recovery system are formidable. Studies into the control system requirements for 
these systems are just beginning [10][12]. 

The final consideration in determining the power requirements is in the choice of rf 

frequency. In general, it is preferred to have as low a frequency as possible without the size 

and capital cost of the accelerator structure becoming excessive. Lower frequency cavities are 

preferred because [13] [22] [24]: 

1. The highest average power rf sources are at 500 MHz and below. 

2. For CW applications, which would become necessary for high-average power electron 
beams, cooling considerations for RT cavities favor lower frequencies. Although the heat 
losses are larger at lower frequencies, the larger cell surface area provides for more 
effective heat transfer. 

3. The larger surface area per cell allows for more effective and simpler cooling of the cell 
walls. 

4. Lower frequency cavities have the potential to produce higher quality electron beams for 
three reasons. First, for a given spatial length of each electron micropulse, it will span a 
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smaller A(J» phase of the rf accelerating fields which leads to more uniform acceleration 
properties. Second, the transverse dimensions (apertures) of the cavity are larger so the 
"transverse wake-field effeas," or beam emittance induced in traveling through an 
aperture due to the image charges formed on the surface, is reduced. Third, the cells are 
physically longer since the length of a cell is approximately Therefore, for a given 

acceleration gradient, the beam will pass through fewer apertures to reach a given energy. 

5. The volume of the cavity is proportional to so a reduced frequency leads to an 
increased cavity volume and higher stored energies. The electron beam traveUng through 
the cavity creates higher order modes (HOM) by virtue of its charge. In higher energy 
density cavities, these HOM will disturb the cavity fields less and reduce the threshold of 
beam-breakup instabilities that can result. 

6. Individual cavities can be more easily electromagnetically uncoupled which allows for 
better stabilization of the fields within. 

The disadvantages of lower frequency operation are the increased size of the strucmre, 
and the fact that acceleration gradients must be lower. They are also more difficult to 
manufacture because of the increased size. 

Figure 12 is a plot of equation 3.6 given values of = 10 MV/m, W = E x = 96 
MeV, Qres = 8x10^, and T = 4.2 K. The minimum occurs at a frequency of 300 MHz. If 
or T decreases slightly, then this minimum will shift to higher frequencies. Therefore, a 
frequency of 400 MHz is a suitable design point for the SC accelerator. The wall losses for the 
SC structures will be approximately = 1.2 kW. This value is insignificant to the overall 
power efficiency of the system. However, this value wiU determine the requirements for the 
liquid helium refrigeration support system. 

The RT accelerator operating frequency is not as easy to determine. There are three 
criteria which determine the possible operating frequency range. First, the maximum field 
strength, £max- allowed within a cavity before sparking occurs is £„ax ~ 1.5£/c/>/DA/ [23] 
where DM is an appropriate design margin and E^^ (MV/m) is the Kilpatfick sparking criteria 
field strength determined from the trancendental relation 

/ = 1.64 £^P (3.9) 

A design margin of four is generally used. 

Second, a minimum field strength within the cavity is necessary to ensure adequate 
energy storage so that the energy removed in accelerating an electron pulse is a small fraction 



35 



sc - Wall Losses 




Figure 12 - SC accelerator wall losses 



of the total cavity stored energy. This minimum field strength, is given by [23] 

_ ^ I 2rt/ M f^yg 
" 77 (AUlU) 



(3.10) 



where AUlU is the fraction of energy removed from the cavity by an accelerating electron 
pulse and M represents the pulse repetition rate in integer numbers of rf cycles. A value of 
M = 12 means that an electron pulse is accelerated every 12'^ rf cycle. 

The final criteria needed to determine the operating frequency for RT structures is the 
maximum allowed wall losses, Equation 3.7 establishes the relationship between and 

/ for a fixed . 

The relationship between these three criteria is shown in Figure 13(a) with values 
k^ = 3.5 = 8xl0'^ M = 12, = 500 mA, AU jU = 0.02, and 

^max = ^ curve Can be lowered only by changing M given the design 

constraints of our problem. Current accelerators operate with values of 12 < M <20. For a 
given value of M and AUlU , the wall losses are less at lower frequencies because of the 
reduction in field strength. A value of A/ = 12 was chosen to be consistent with current 
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Acceleration Gradient Limitations 



RF - Wall Losses 




(a) (b) 

Figure 13 - RT accelerator (a) field strength requirements, and (b) wait losses 

technologies. As M is reduced, the peak currents are reduced as is the excitation of HOM 
within the accelerator cavities. However, less time is available to damp these modes out. 

Within the design constraints of the PEL. a maximum rf frequency is desired to limit tJie 
wall losses in the RT structure. A maximum acceleration gradient is also desired in order to 
reduce the active length of the accelerator structure. However, this is limited by the resultant 
increased power losses. Therefore, an operating frequency of 200 MHz is chosen for the RT 
accelerator. From Figure 13(b). the wall losses are Pi = 2.9 MW with an £<, = 1.5 MV/m. 

This analysis does not address the suitability of a structure with these parameters for use 
in an FEL. Extensive modeling and numerical simulations would be required to determine the 
beam quality and beam transport characteristics that could be achieved with such a structure. 
However, these parameters have been validated by a comparison to an actual 180 MHz cavity 
design [12]. 

It is prudent to note that the trend for future compact FEL accelerators is to include 
recirculation with energy recovery where a beam is bent around to pass through the accelerator 
two or more times to reach its final energy (2][13][25j. This allows accelerators to be more 
compact in that the active length of the accelerator is smaller by a factor approximately equal 
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to the number of passes. The power estimates for these specific architectures will be nearly 
identical to the analysis presented here. However, the control system issues are more complex. 

C. ELECTRON INJECTOR 

The purpose of the electron injector is to produce high density pulses of electron current 
and inject them into the accelerator at the proper phase relation to the accelerating rf fields. 
The injector consists of two major components: the electron gun and the buncher. The buncher 
systems contain rf cavities similar to accelerator cavities. These cavities are designed to 
produce bunched electron pulses of low emittance. One of the most promising electron gun 
technologies for the production of high current, low emittance electron beams center on the dc 
photocathode [12] [14], 

Figure 14 shows the functional component interconnections of an injector driven by a dc 
photocathode electron gun. When illuminated by laser light, electrons are liberated at the 
photocathode through the photoelectric effect. A dc high-voltage rapidly accelerates the 
electrons away from the photocathode and directs them into the buncher. The dc voltage 
should be large to prevent the Coulombic space charge effects from inducing excessive initial 
emittance. However, since the electron velocity cannot exceed the speed of light, the voltage 
must be low enough to limit the electron velocity sufficiently below the speed of light when 
the electrons enter the buncher so it can effectively bunch the micropulses before they enter the 
accelerator [22]. A suitable value for the dc voltage is 500 kV. The electrons have been 
accelerated to a kinetic energy of about 500 keV ( velocity, v = 0.85c ) when they enter the 
buncher. 

In describing the electron gun power requirements, it is more convenient to start at the 
high voltage dc electrodes and work back towards the drive laser. The high voltage has a 
practical limit due to dielectric breakdown effects and the emittance limiting effects previously 
mentioned. As buncher and high voltage power supply technology has improved, however, 
this voltage limit has also risen. Currently, and for the near future, these devices can be 
expected to operate at about 500 kV [2][14]. Therefore, for a 500 mA average electron 
current, 250 kW of high voltage dc electrical power must be supplied by the dc electric field to 
the electron current. 

To generate the 500mA of average current, we use a photocathode which emits electrons 
via the photoelectric effect. It is assumed that our accel^ator will be operating in the CW 
mode so that the peak currents in each micropulse is minimized. However, to minimize the 
beam emittance, the electron micropulses should have a spatial length corresponding to only a 
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(a) typical dc photocaihode electron gun 




Accelerator 



(b) electron injector 

Figure 14 - typical electron injection system 



few degrees of rf phase. Reference [22] defines a good injector as supplying a phase spread of 
A(|) S. 10° and a very good injector as having a A<J> ^ 5°. Low emittance electron beams are 

being produced with micropulse phase spreads of approximately A<|> = 5-6° [2][14][26]. For a 
CW 400 MHz system, a A(() = 6° corresponds to a pulse length of 42 ps with an average 
micropulse current of 30 A per rf cycle. Therefore, we require an average electron charge of 
1.25 nC per rf cycle or 1.25xA/ nC per micropulse. At a rf frequency of 200 MHz. we 
require an average electron charge of 2.5xM nC per micropulse. Electron pulses of 3-4 nC. 
peak currents of several hundred amps, and average currents of tens of mA have been 
effectively produced [26] [2] [1 4] with rapid advancements being made. The expectation of a 
500 mA average current in the future is reasonable. 

The photoelectric effect is described by 

=//v + H^^, (3.11) 
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where Ej^ is the kinetic energy of the liberated electron, /jv is the energy of the absorbed 
photon, and Wpg is the work function of the photoemissive materiaL The values of the work 
functions for the different elements range from 3.14 eV (cesium) to about 6 eV [27]. This 
minimum photon energy required to generate a free electron means that photocathodes are 
insensitive to infrared radiation. The amount of charge liberated by a pulse of total energy Ep 
with photons of energy /i v is 



< 7 (nC) = 10 



Epi\iJ)QEi%) 

/jv(eV) 



(3.12) 



where QE is the effective quantum efficiency which is the ratio of the number of liberated 
electrons to the number of incident photons. Table 4 lists the effective quantum efficiency for 
a number of photocathodes under development. 



photocathode 


QE 

(%) 


1 

'^incident 

(nm) 


Comments 


CsKjSb [26] 


0.1 


527 


1-3 nC 


Cu [UCLA] 


.01 


266 


> 1 nC 


Ag [29] 


.03 


248 




Cu+10%BaO [Aganofov] 


.17 


266 




LaBe[31][32] 


.06 


349 


T = 1500 K 


GaAs [14] [2] 


18 


527 




CsTe [33] [34] [35] 


13 


263 


long lifetime 


CsTe [34] [35] 


19 


251 





Table 4 - Photocathode quanmm efficiencies 



It is reasonable to expect that a photocathode will be developed suitable for a shipboard PEL 
with a quanmm efficiency of QE = 20%. However, a QE = 5% would be a more 
conservative expectation, and would account for QE variations over the lifetime of the 
photocathode. From 3.12 we can show that the required laser input power to the photocathode 
to generate a required current is 
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(W) 



(3.13) 



P = 



/tv(eV) 
10 QE{%) 



/(mA) 



For a frequency doubled Nd;YAG laser (X ~ 532 run) and a QE ~ 5%, an average current of 
500 mA would require an input power of 1 1.7 W incident on the photocathode. 

The frequency doubler is a non-linear electrooptic crystal in which second harmonic 
generation occurs. If we align the crystal’s crystallographic axes such that the index of 
refraction for both the fundamental and second harmonic frequencies are equal, the crystal is 
"phase-matched". Under the phase-matched condition, two photons of the fundamental 
essentially combine to form a new photon of twice the frequency, that is 

— > k 2f 



H(£> + luo = H(2a>) 



(3.14) 



The efficiency of frequency doubling, ranges from 20-65% depending on the input power level. 
Present photocathode electron guns operate with a frequency doubling efficiency, of about 
50% [26][14], but that can vary depending on the actual incident power intensity. For linac 
applications, the length of the electro-optic crystals are limited by group velocity mismatch 
which tends to increase the pulse length and limit the frequency doubling efficiency. Using 
type 1 phase-matching, however, efficiencies of more than 60% can be obtained [36]. Due to 
our high-power requirements, a frequency conversion efficiency of T[fj = 55% is assumed. 

An optical amplifier is a pumped laser medium which amplifies the input signal through 
stimulated emission. An optical amplifier is described by the expression [37] 

ln(G,^) + -^[G,p-1] = (3.15) 

where is the unsaturated gain coefficient, d is the single pass length of the amplifier 
medium, G^p is the overall single pass gain, is the saturation photon flux density, and <1H0) 

is the input photon flux density. is a function of only the amplifier medium while is a 
function of the amplifier medium and is directly proportional to the population inversion per 
unit volume. Therefore, Po wiU determine the power requirements for the amplifier. 

To maximize the power efficiency of the amplifier, the amphfying medium must be 
efficiently pumped and multiple passes must be made by the light wave through the amplifying 
medium. To efficiently pump the medium, laser diode pumping provides the most efficient 
method because the pump energy is concentrated at the pump bands of the Nd lasing medium. 
The width of the diode pump spectrums are on the order of 5 nm FWTHM [38]. Furthermore, 
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diode lasers, driven by a low dc voltage, can achieve a quantum efficiency of near unity 
[39][40]. One particular diode-pumped optical amplifier design utilizing a confocal, parabolic 
mirror cavity with an internal Brewster plate to increase the path length of the beam within the 
crystal, has achieved small signal gains of 38 dB with a pump to optical power efficiency of 
40% for the Nd:YLF medium and 25% in the Nd:YAG medium operating in the CW mode 
[41], Amplified spontaneous emissions are only about 1% of the total output power. The 
Nd:YAG amplifier, however, was able to support much more pump power than the YLP 
amplifier which fractured at lower pump powers. For the estimation of amplifier power 
requirements, it is assumed that an optical amplifier for my FEL photoinjector has a pump 
power amplification efficiency, , of 25%. 

Nd:YLF and Nd;YAG mode-locked CW lasers are commonly used as photoinjector drive 
lasers because they are relatively rugged, produce short pulses (about 75 ps for YAG [39] and 
50 ps for YLF [42]), have large optical gain per unit distance, and can be easUy amplified. 
The large optical gain means they have a low threshold power required for CW operation. 
With laser diode pumping, we can achieve much higher pumping efficiencies. The resulting 
inefficiencies are predominantly due to incoupling losses, the fraction of energy absorbed by 
the lasing medium, and the spontaneous emission losses. Optical-to-optical efficiencies of the 
diode-pumped laser have been demonstrated at 25% for Nd lasers in CW mode [43]. Diode 
lasers delivering 60 W for a 200 psec duration are commercially available. Advances in the 
development of higher-powered diode lasers and improvements in the coupling optics are 
expected to lead to pulse durations below 10 nsec in the near future [43]. This is sufficient to 
support the 42 ps pulse lengths needed to support this FEL. 

Current Nd drive lasers used in photocathode electron guns have self phase modulation 
(SPM) pulse compression systems with an incoupling loss of 25% [42]. These guns drive 
accelerators operating at rf frequencies of several GHz where narrow electron pulses are 
required. The pulse compression systems do not require input biasing or pumping power but 
do limit the peak output power from the laser to about 2 W. Pulse compression systems, 
however, should not be required for the low accelerator frequencies under consideration for a 
high-power FEL. If we assume a laser ouqjut power of 2 W at a waU plug efficiency of 25%, 
only a few watts of power are needed. Therefore, the drive laser power and pulse compression 
losses can be neglected from the prime power analysis. 

The power requirements of the electron gun now can be summarized in terms of Figure 
15. The total required low voltage dc power, Puc< is given by 
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Low Voftao* cJc 600 kV dc 

170 W 260 kW 




Figure 15 - Power flow in the dc photocathode 



^Ldc 




1 /tv(eV) 
10 < 2 £(%) 






(W) 



(3.16) 



and the total high voltage dc required. P/y*. is the average current times the voltage applied. 
For a 500 mA average current and 500 kV acceleration, we will require Pt^c = 170 W and 

= 250 kW. 

The buncher cavities are RT cavities whose geometry and fields are optimized to provide 
bunching of the electrons into micropulses with low emittance and energy spread. SC cavities 
and their associated geometries are not efficient for tliis purpose so RT cavities are used 
exclusively. The electrons liberated from the photocathode are bunched longitudinally and 
focused to a waist as they are injected into the accelerator. Solenoidal magnetic fields are also 
used to assist the bunching process. 

For the defined architecture and parameters, the buncher cavities must provide 3 MeV 
acceleration to the electrons during this bunching process. Since the electrons are entering tlie 
buncher at a velocity of v = 0.85c. the fields must be rather su^ong to bunch the elecu^ons. To 
determine the wall losses in the buncher, Punj< the same model used for estimating the 
accelerator wall losses is employed. To achieve 3 MeV acceleration, three cavity cells at a 
frequency of 200 MHz and Eg = 1.33 MV/m is assumed. To account for added losses due to 
the required geometry changes, the shunt resistance constant is assumed to be A , = 1.9 
(M£2/m/MHz Therefore, we can estimate ~ 150 kW. At a frequency of 400 MHz. 
only 100 kW of power would be dissipated. 

Although electromagnets are often used to produce the solenoidal fields, large permanent 
magnets with small tuning electromagnets would be more appropriate for high-current 
applications. Therefore, they will not add significantly to the power requirements. 
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D. RF SYSTEMS 



Electron accelerators have historically operated at rf frequencies of more than 1 GHz, low 
average power, and high peak powers in a pulse power electron beam format. Klystrons have 
been the only rf sources suitable for these operating conditions. For rf frequencies of 200-500 
MHz, however, additional rf sources may be practical. 

Table 5 is a listing of several of the most powerful CW klystrons currently on the 
market. 





TH2089 


YK1350 


YK1303 


K3513 


Units 


Reference 


[44] 


[44] 


[44] 


[45] 




frequency 


352 


352 


508 


352 


MHz 


input 

voltage 


87 


90 


90 


100 


kV 


input 

current 


17.1 


16.8 


18.2 


20 


A 


Power 


1 


1 


1 


1.3 


MW 


efficiency 


68 


66 


61 


65 


% 


length 


4.75 


3.75 


3.75 


4.2 


m 



Table 5 - Commercially available super-power CW klystrons 



Klystrons have evolved over several decades so revolutionary advances in their power 
capabilities within the next 20 years are not anticipated. Some of these tubes are designed to 
operate in a vertical orientation while others are designed to operate in a horizontal orientation. 
High-power klystrons are available in either configuration. The output power and power 
conversion efficiency of a klystron for a shipboard FEL can be expected to have an output 
power of about 1.5 MW at an efficiency of 65%. 

Some alternative rf sources for accelerator applications are listed in Table 6. The high- 
power entries in this table are projections and have yet to be developed. However, these 
alternate technologies promise rf sources that operate at lower voltages, higher efficiency, and 
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MBK 


MBK* 


klystrode 


klystrode* 


lOT 


Units 


Reference 


[46] 


[46] 


[47] 


[48] 


[49] 




frequency 


425 


500 


267 


476 


476 


MHz 


input 

voltage 


17 


39 


65 


25 


66 


kV 


input 

current 


11 


41 


5.5 


70 


9.3 


A 


Power 


0.1 


1 


0.25 


1.3 


0.4 


MW 


efficiency 


55 


64 


70 


75 


> 65 


% 


length 


1.35 


2.5 


NA 


NA 


2.1 


m 




400 


800 


NA 


NA 


500 





Table 6 - Alternative CW rf power sources (* - denotes projection and is not 
commercially available) 

have a reduced component weight and volume. It is expected that the multiple beam klystron 
(MBK) efficiency can be increased to >80% with further development [24], The klystrode and 
inductive output tube (lOT) are essentially the same device with klystrode being Varian 
Associates trademark for this device. It has the disadvantage of a smaller gain. Therefore, a 
larger rf drive is required for the device (several kW for a megawatt device). This will add to 
the complexity of the overall rf system. 

E. SUPPORT SYSTEMS 

A wide variety of support systems will be required for a shipboard PEL. Feedback 
control systems will be required to maintain rf power source phase and amplitude stability. The 
2.5 MW of heat energy dissipated in the beam dump and the 3 MW of heat energy dissipated 
in the RT structure as waU losses will need an appropriate cooling system. System monitoring 
and fire-control interconnections will also place a load on the ship’s electrical system. These 
systems, however, can easily be adapted to existing shipboard systems. Current steam 
propulsion plants dump heat energy well in excess of 50 MW at a cost of tens of kW of 
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seawater pump power. Therefore, a reasonable estimate is that these systems will add an 
additional 100 kW to the ship’s electrical load. 



There are three support systems, however, that need to be addressed separately. First, a 
SC structure will require a liquid helium refrigeration system. The power requirements for this 
system can be estimated by [17] 



where Pn is ship’s service AC power consumed by the refrigeration system, Ta is the ambient 
(room) temperature (7^ » 300/f), Tne is the liquid helium boiling point temperature 
(7^^ ~ 4.2^), and is the refrigeration efficiency. The second term in 3.17 represents the 
Carnot thermodynamic efficiency of the system and the third term is the refrigeration system 
(compressor, etc.) efficiency. Typical values of t\k are approximately 25-35 percent [14][16]. 
For a 30% efficient refrigeration system removing a 7^ = 1.2 kW, our shipboard FEL will 
require 7^ = (235x7^ ) » 250 kW of refrigeration power if the laser were to operate 
continuously. This is not only a significant amount of power but the size of the refrigeration 
plant would be large [2] [14]. 

The refrigeration plant for this shipboard FEL, however, can be considerably smaller 
since it will not operate continuously. The rate of liquid helium boUed off during FEL 
operation is 



where V is the rate of helium boiled of during laser operation (gal/sec), A is the atomic weight 
of helium (4 g/mole), p is the density of helium at Tne (7.62 Ib/ft^), and is the latent heat 
of vaporization of helium (84 J/mole). The operational requirements specify that the FEL must 
operate for 150 seconds over a 20 minute interval. Conservatively assuming that the laser 
operates continuously for 150 seconds with a wall loss of 7^ = 1.2 kW, 18.5 gal of liquid 
helium is boiled off. This is a relatively small amount of required stored liquid helium. 
Therefore, the refrigeration system can be designed to recompress and refrigerate this boil-off 
over a longer period. 

The refrigeration system sizing criteria, therefore, is based on two factors. First, the 
refrigeration system must be able to remove the latent heat of vaporization from the helium 
boiled off during laser operation within a specified time period. Arbitrarily choosing this time 
period as 2 hours, then the refrigeration power required is 




(3.17) 



V = Pl 
P <?v 



(3.18) 
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Pyj = 250 kW (100 sec/2 hrs) = 3.5 kW. Most importantly, the refrigeration system size and 
weight has been reduced significantly. A reserve storage capacity of liquid helium can reduce 
this requirement even further. The second factor is the time required to reduce the temperature 
of the accelerator from room temperature to operating temperature following an import 
maintenance period. To perform this cooldown within 48 hours would require a refrigeration 
system that would be excessively large for the shipboard environment [14]. A shore facility 
could perform this cooldown prior to a warship getting underway and eliminate the need for 
this sizing criteria. 

The second support system that needs to be addressed is the network of electron beam 
bending magnets which control the beam trajectory. Large permanent magnets would also be 
appropriate for this purpose. Although they are much more expensive, they require no external 
electrical power source. If electromagnets were used and a magnet failed, the beam would 
impact the wall of the beam tube creating a large radiation hazard for adjacent personnel and 
causing a shutdown of the PEL. 

The final support system deals with the optical resonator. The resonator will have 
megawatts of power circulating between the mirrors. Special mirror configtiradons have been 
designed for high-power applications [2] to reduce the power intensity on the mirrors and the 
onset of damage. The optics will need, however, an active stabilization system to maintain 
alignment during ship vibrations, flexing, and torsional movements. Furthermore, neglecting 
outcoupling of the laser light, the mirrors will also absorb some energy and need a cooling 
system. Mirrors have been developed with a reflectance of 99.975% at X = .87 |im [50]. In 
general, mirror reflectivity improves as the optical wavelength moves into the infrared. 
Assuming a reflectance of 99.9% with a resonator stored power of 5 MW, 5 kW of energy is 
absorbed in the mirrors. Suitable stabilization and cooling systems have already been developed 
[24][51]. The power requirements can be included with the overall support system power 
allocation, of 100 kW. 

F. SUMMARY OF POWER REQUIREMENTS 

The power requirements for the PEL architecture of Figure 10 under the design 
constraints of a naval application are summarized in Table 7. Since about 95% of the electron 
beam power is recovered in this architecture, the listed power requirements are nearly 
independent of the output power level. For example, the accelerator wall losses are solely a 
function of the stored cavity field strength. Therefore, the wall-plug efficiency of the FEL 
scales as the output power. The rf source efficiency is based on commercially available 
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parameter 


SC value 


RT value 


units 


/ 


400 


200 


MHz 




100 


100 


MeV 


lavg 


500 


500 


mA 


Ea 


10 


1.5 


MV/m 


z 


311 


990 


Q/m 


Qres 


X 

o 






Qo 


3x10’ 


5x10" 




E iMc 


350 


350 


W 


PhcJc @ 500 kV 


250 


250 


kW 


Pinj 


2 


2 


MW 


E Unj 


150 


150 


kW 


P 

‘ beam 


50 


50 


MW 


Tlo 


2 


2 


% 


Pout 


1 


1 


MW 


Eret 


49 


49 


MW 


Arec 


95 


95 


% 


p 

‘ dump 


2.45 


2.45 


MW 


Prec @ ^rec = -95 


46.6 


46.6 


MW 


Pl 


.0012 


2.9 


MW 


Pr @ Or = -3 


275 




kW 


P RF @ Orf = -65 


6.6 


9.2 


MW 


Pss 


100 


100 


kW 


wall-plug efficiency 


17.4 


10.4 


% 



Table 7 - Summary of the power requirements for a shipboard PEL 



klystrons. Alternative rf sources promise an increased efficiency but require further 
development before they can produce high-average power. If extrapolations are made to higher 
power for these devices, then the PEL total power requirements can be reduced by about 10%. 

Por the RT accelerator, it was necessary to define the pulse repetition interval represented 
by the parameter M to determine the accelerator operating frequency and wall losses. Prom 
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equation 3.10 and Figure 13 it can be seen that reducing the value of M can reduce the RT 
cavity wall losses. However, the optical resonator must have dimensions set so that when an 
electron pulse passes through the undulator, it overlaps with an optical pulse and can transfer 
energy to the optical field. Therefore, the mirrors must be separated by a distance 
L = nMcIlf , where n is an integer and equals the number of optical pulses stored within the 
resonator cavity. For the values in Table 7 and assuming A/ = 12, L = 9n meters for the RT 
accelerator parameters. If the resonator cavity was 18 meters long, then there would be two 
optical pulses circulating within the resonator. 

For the SC accelerator based FEL, M is not an important parameter in determining power 
losses or the rf frequency because of the high cavity Q and large accelerating field strengths. 
However, the pulse repetition rate will affect the length of the resonator in the same way as for 
the RT accelerator based FEL. For the SC parameters defined in Table 7, and assuming 
M = 12, the resonator length is L = 4.5n meters. 
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IV. PRIME POWER DISTRIBUTION 



In the previous chapter, the power requirements for a megawatt FEL laser weapon for 
shipboard applications were outlined. The two components that burden the prime power 
distribution system are the dc photocathode high voltage and the rf sources for the electron rf 
accelerator. It was determined that the dc photocathode would require 250 kW at 500 kV dc 
while the rf sources would require 6-10 MW at 80-100 kV dc. This amount of power caimot 
be delivered through the conventional ship’s service power distribution system. In this chapter, 
the characteristics of a shipboard prime power distribution system capable of delivering power 
to these FEL components are explored. 

The prime power supply requirements are summarized in Table 8. These requirements 
are established to be consistent with the FEL power requirements of Chapter III. 





rf sources 


dc photocathode 


voltage 


100 kV 


500 kV 


power 


10 MW 


250 kW 


type 


dc 


dc 


duration 


15 sec 


15 sec 



Table 8 - Prime power requirements to support a high-power FEL 

The value of 100 kV for the rf source input voltage was chosen to support the high-power 
klystron operation. This is a limiting value and would likely be about 40 kV when alternate rf 
sources are developed. 

A. POWER GENERATION AND DISTRIBUTION 

The FEL system wiU require a dedicated distribution bus. Therefore, this bus is not 
limited to the normal ship’s service distribution system operating parameters of 450 V, 3 
phase, and 60 Hz. Several studies have been performed to outline the prime power system for 
high-power, electricaUy-powered weapons with power requirements similar to FELs 
[52][53][54][55][56], The power generation and distribution aspeas of these studies are equally 
pertinent to this discussion and are summarized in this section for continuity purposes. 
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Follow-on sections will address components specific to the FEL systems. 

The prime power distribution system is illustrated in Figure 16. Power generation and 
distribution component weights and volume are summarized in Table 9. Power is diverted from 
a propulsion gas turbine engine to an auxiliary generator to power the FEL distribution bus. 
Modifications to the pinion that connects the turbine to the main reduction gears have been 
designed so that the turbine power is automatically diverted to the auxiliary generator without 
operator action when load is drawn by the FEL loads [52][53]. Simulation smdies have shown 
that this diversion of power will only cause a momentary reduction of ship speed of a few 
percent when operating at flank speed over the duration of the laser shot [53]. The lED 
generator pictured in Figure 16 represents the electrical generator which would power an 
electric drive ship. This element would be deleted for a ship powered with mechanical drive 
where the propeller shaft is directly coupled to the reduction gears. 

The FEL distribution bus would operate at about 5 kV, 3 phase, 800 Hz. As bus 
frequency is increased, the size and weight of generator and transformer components are 
reduced. However, their efficiency is also reduced. A frequency of 800 Hz was determined to 
be the best compromise between these two competing factors [52][55]. The 5 kV bus voltage 
was chosen due to conductor sizing limitations. A low bus voltage causes the transmitted 
current to be high resulting in the need for large cable conductors. Similarly, a large voltage 
results in increased cable size due to insulation requirements. 

B. 20 MVA TRANSFORMER 

Computer algorithms have been developed for the prediction of transformer weights for a 
pulsed, high-power electrical load. Over a power range of 10-50 MW, the weight of an 
adiabatic (no forced cooling) transformer is given by the expression [57] 

Weight = 0.0505 * [0.693 + 0.307 (F/25)“°-^^] 

* [0.931 -I- 0.069 (Fo/100)'-^] * [0.242 -i- 0.758 /‘^•^^^] (Ibs/kW). (4.1) 
For a 1-5 MW transformer, the weight is given by 

Weight = 0.1275 (r/120)®-^‘ * [0.693 -i- 0.307 (F/25)^’^] 

* [0.931 -I- 0.069 * [0.242 -i- 0.758 (ibs/kW) (4.2) 

where r is the total run time (sec), P is the power level (MW), is the dc load voltage (kV) 
assuming rectification without phase-control at the output of the transformer, and / is the bus 
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Figure 16 - Prime power disU'ibution system. After Ref. [52] 
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component 


weight 

(kg) 


volume 

(m^) 


gas turbine 


20,400* 


46* 


Main reduction gears 


6800* 


3.8* 


auxiliary generator 


8000 


4.3 


transmission cable 


4500 


1.0 


auxiliary support 


3300 


NA 



Table 9 - Summary of power generation and distribution component weights and 
volume. * - part of propulsion system and does not add to PEL system 
totals; NA - not available 

frequency (kHz). These equations assume an input voltage of 1-5 kV and have been validated 
for output voltages of 20-200 kV, run times of 1-120 seconds, off times of 2-300 seconds, and 
frequencies up to 2 kHz. The studies which produced these relations also showed that adiabatic 
transformers provide for the lowest size and weight transformer charaaeristics for pulsed- 
power applications. 

In order to conform with Navy standards, all components must be derated. The derating 
factors for both current and voltage are 0.7. The transformer which will supply the 10 MW 
bus for the rf sources was derated by (0.7)". Therefore, it was rated at 20 MW. From equation 
4.1, the required 100 kV, 20 MW transformer would weigh 620 lbs (280 kg). The weight of a 
100 kV, 2 MW transformer would be 157 lbs (71 kg). Therefore, there is a definite weight 
advantage to using a single transformer to supply a single 100 kV bus instead of using a few 
smaller ones. 

A volume determination can be made from the method described in Ref. [55]. Figure 17 
in Ref. [55] gives the volume of a transformer as a function of frequency and power level. 
From this curve, the 20 MW transformer would occupy approximately 5 m^. 

C. 100 KV RECTIFIER 

The rectifiers were assumed to be composed of silicon-controlled reaifiers (SCRs) instead 
of diodes. This is because the accelerator rf sources and dc photocathode voltages need to be 
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controlled to tight tolerances. Use of a phase-controlled rectifier allows for a control system to 
regulate the output voltage. The weight of the 100 kV reaifier was estimated to be 1500 kg 
while the volume was estimated to be 8 m^. The analysis used to determine these values is 
presented in the Appendix. 

D. 500 KV POWER SUPPLY 

The 500 kV power supply must provide an average current of 500 mA. A voltage 
multiplication circuit would provide better performance [12] and be smaller and lighter than an 
equivalent transformer and rectifier combination. The weight and volume of a 500 kV voltage 
multiplier was estimated to be 4100 kg and 3 m^, respectively. The analysis used to determine 
these values is presented in the Appendix. 

E. ENERGY STORAGE 

An important consideration for any weapons system is how it operates under casualty 
conditions. A dedicated energy storage element could provide a limited capability for PEL 
operation should the normal power distribution be interrupted. Several methods of energy 
storage were considered; storage batteries, superconducting magnetic energy storage, flywheels, 
and capacitor banks. Only capacitor banks are adequate to deliver megawatts of power within a 
fraction of a second and yet be compatible with the shipboard environment. The weight and 
volume of high-energy storage capacitors is proportional to the storage capacity. Modem 
capacitors can store 3 KJ/kg [53] and it is anticipated that 4 KJ/kg is achievable with some 
development [54][58]. These capacitors would occupy a volume of 5 MJ/m -^ [53][58]. 

The run time of 15 seconds for the PEL is based on the time available to engage a target 
and not on the laser dwell time required to destroy the target. Therefore, it may not be 
necessary to store enough energy for a full 15 second engagement. 20 MJ of stored energy 
would provide for a limited 2 second engagement. This amount of stored energy would 
require a capacitor bank that weighs 5000 kg and occupies 4 m^. This capacitor bank would 
also require a power electronics package on the order of the 10 MW rectifier system in size 
and weight. The detailed estimates of this package were not made. 

F. SUMMARY OF FEL PRIME POWER DISTRIBUTION 

The weight and volume estimates for the prime power distribution system is summarized 
in Table 10. Those components which are normally part of the propulsion train are not 
included. The transformer and rectifier values were determined for a 100 kV FEL bus voltage. 
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component 


weight 

(kg) 


volume 

(m^) 


power generation 


15800 


5.3* 


transformer 


280 


5 


rectifier 


1500 


8 


500 kV power supply 


4100 


3 


totals 


21680 


21.3 



Table 10 - Prime power weight and volume estimates. * - does not include 
auxiliary support equipment 

These weight and volume estimates can be reduced appropriately as alternate rf sources 
become available. 

Energy storage is not included in this summary. However, its inclusion would result in 
an increase of more than 5000 kg and 5 m ^ to the totals. This would be a significant increase 
in overall system size and weight. 
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V. HIGH-POVVER FEE OPERATION 



In the previous chapters, the design constraints and prime power considerations for a 
shipboard high-power FEL have been outlined. Various electron beam parameters have been 
defined so that an FEL could provide a 1 MW optical output. In this chapter, the conversion 
of electron beam energy to optical energy is explored. The undulator parameters considered in 
this chapter are not meant to specify a required design but to demonstrate that the defined 
electron beam parameters can. in fact, produce the desired output. Both the undulator and 
optical resonator parameters are also constrained by the overall system goals of minimum size 
and weight. 

The electron beam parameters that must be considered to adequately characterize the FEL 
operation are summarized in Table 11. 



parameter 


value 


units 


Beam energy, 


100 


MeV 


Peak Current. Ip^^k 


360 


A 


M 


12 


- 


Micropulse length, (400 MHz) 


42 


psec 


Micropulse length, (200 MHz) 


84 


psec 


Beam radius, 


330 


pm 


Energy spread, Ayly 


0.002 


- 


e (rms) 


0.1 


mm-mrad 



Table 1 1 - Electron beam parameters 

The electron beam radius, the electron beam energy spread, Ayly, and the electron beam 
emittance. e. are estimates made from existing accelerator electron beams. A high quality 
electron beam from a 1.5 GHz accelerator with a peak current of 200 A has the values 
r* = 164 p.m. Ay/y = 0.1%, and an e = 0.1 mm-mrad [2]. The electron beam radius and 
emittance scale approximately as the current [12]. It has already been established that emittance 
is reduced in lower rf frequency accelerators. Furthermore, electron energy spread Ay is 
conserved as an electron beam is accelerated [22]. Therefore, values of r^, = 330pm, 
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Ay/y = 0.002, and e = 0.1 nun-mrad are reasonable assumptions. 



A. FEL OPERATION ■ 3 pm WAVELENGTH 

The operating wavelength of an FEL is determined by the resonance condition, equation 
2.6, where K - eBgXollninc^. For operation at A. = 3 pm, a range of FEL undulator 
parameters are summarized in Table 12. These parameters were defined in Chapter II, sections 
A and D. 



(cm) 


K 






Bg (peak) 
(kG) 


1 


4.7 


.48 


3.3 


71 


3 


2.6 


.43 


1.9 


13 


5 


1.9 


.39 


1.5 


5.7 


10 


1.1 


.28 


.96 


1.7 



Table 12 - Undulator parameters for FEL operation at A. = 3 p;/j 



The undulator field strength is limited by technology. The maximum peak undulator field 
strength that can be produced is given by the expression [59] 

( 5.47 - 1.8 ) 

B„=33.3e “ (AG) (5.1) 

where g is the gap between the pole pieces of the undulator magnets. The strongest field 
strength that has been manufactured in an undulator was Bg = 18 kG with a glXg = 0.1 1 [60]. 
The minimum gap distance is determined by the transverse dimension of the optical mode. 
Under these limitations the minimum undulator wavelength that can be considered for this high 
power FEL is « 3 cm with a g ~ 5.4 mm. 

The electron micropulse passing through the undulator is several hundred slippage 
distances long. Hence, the optical pulse will be approximately the same length as the 
micropulse and the duty cycle between pulses will also be equal. For the parameters of Table 
11, the duty cycle, d, will be about 1/720. The optical pulse amplitude and phase 
characteristics will be periodic along the length of the pulse with a period of about one 
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slippage distance. Therefore, simulation of the FEL electron-optical interaction over one 
slippage distance is sufficient to model the entire pulse. 



The relationship between average power output of the optical resonator and the complex 
field strength is given by 



r> 2 jC 

Pavg - d — 



y~mc" 



AneNlK^i^) 



\a\^ (kW) 



(5.2) 



where Q is the quality factor of the optical resonator. Over the length of the undulator, 
however, the complex field strength and optical mode cross-sectional area are not constant due 
to diffraction effects. This variation is a direct function of the Rayleigh length, For a 
high-power FEL, we want and vt’^ to be as small as possible to minimize the power density 
on the resonator optics. However, if < r ^, , then the gain of the FEL is reduced since not aU 
of the electrons can interact with the optical mode. As explained in Chapter 111, section E, this 
is accounted for by setting j = j F , where F = when simulating the wave equation. 

The value h„, however, is dependent on the Rayleigh length of the optical resonator. 
Assuming Gaussian beam propagation, the value Zg ~ LI^Yt. minimizes the volume for the 
optical mode over the length of the undulator and provides for maximum gain when Wg > rf, 
[5]. With = L/VI2, we can use equation 5.2 and neglect the variation in lo I. 

The optical parameters for several undulators are listed in Table 13. The Rayleigh length 
of the resonator is assumed to be Zg = L/Vl2. 



(cm) 


N 


L 

(m) 


(mm) 


j 


3 


15 


.45 


.35 


6.3 




25 


.75 


.46 


17.6 


4 


15 


.6 


.41 


6.3 




25 


1 


.53 


17.4 


5 


15 


.75 


.46 


6.2 




25 


1.25 


.59 


17.1 



Table 13 - Optical parameters for FEL operation at = 3 (X/// 
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Several simulations were performed to determine the relationship between j , \a I, and Q. The 
results of these simulations for values of j = 4, 7, and 13 are shown as curves in Figure 17. 
Overlayed with these curves is the la I vs. Q relation determined by equation 5.2 for a 1 MW 
optical output and an undulator with = 3 cm, N = 16. and j = 7. 



required la| lor given Pavg output 




Figure 17 - la I vs. Q dependence for several values of j. The curve marked 
"P = 1 MW" is the la I vs. Q relationship for an undulator with j = 7, = 3cm, 

and = 16. The intersection of the two solid curves is the required operating 
point 



The optical mode within the resonator must have an la I = 17 and the resonator must have an 
outcoupling of 20% (Q ==4.5). From Figure 17, it is obvious that more than 1 MW can be 
exuacted from the elecu-on beam by increasing the resonator Q. However, this would result in 
an undulator extraction efficiency of more than 2%. 

The average power within the resonator cavity is 5 MW and is contained within an 
optical mode size of = 0.364 mm. Assuming a Gaussian profile, essentially all of the 
power (£ 100 W) is contained within a cross-section of 8tv„ = 2.7 mm at the ends of tlie 
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undulator where the spot size is largest. Allowing for alignment tolerances, the undulator gap 
distance g = 5.4 mm is adequate. 

B. ANALYSIS OF FEL OPERATION 

The FEL saturates at an average power of 1 MW for y = 7 and Q ~ 4.5. The steady 
state characteristics of the FEL interaction are shown in Figures 17, 18, and 19. Figure 18 is 
the gain and phase curve, as defined in section H of Chapter II, for this FEL near saturation. 
The undulator’s sensitivity to electron beam quality, as explained in section G of Chapter II, is 
proportional to the number of periods. Therefore, the gain curve is virtually insensitive to 
beam quality since the undulator only has 16 periods. Figure 19 is the phase space evolution of 
an electron pulse at saturation. Approximately 60% of the electrons are trapped in closed 
orbits. The trapped electrons leave the undulator with phase velocities concentrated about v ~ 
-5 while the non-trapped electrons exit with phase velocities spread about v = 5. 




- 0.27 



0.27 



0.0 
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Figure 18 - Gain spectrum at saturation for operation at A. = 3 pm 
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-rt/2 ^ 3rt/2 Ox 1 



Figure 19 - Phase space evolution at saturation for operation at \ = 3 pm 

The startup characteristics of the FEL are illustrated in Figure 20. The evolution of the 
optical mode over n = 150 micropulses beginning from spontaneous emission is shown. The 
left window shows the complex field strength over one slippage distance, tlie center window 
shows the power spectrum, and the right window shows the electron distribution. The electron 
distribution of this simulation agrees with the phase space distribution of Figure 19 in that the 
electrons are concentrated in two distinct phase velocity groups. The power spectrum has a 
single mode and so there is no modulation of the complex field strength envelope. The 
steady-state complex field strength is l« I = VP (n ) « 17. The optical field reaches saturation in 
about 50 passes (1.5 psec for an rf frequency of 400 MHz or 3 psec for an rf frequency of 
200 MHz) for the FEL to reach the steady-state power level of 1 MW. In a combat 
environment, this is essentially an instantaneous engagement. 
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Figure 20 - Optical field characteristics at saturation for operation at A, = 3 jam 

C. FEL OPERATION - 10 |ini WAVELENGTH 

In the maritime environment, there are two principal wavelength bands in the infrared 
(IR) suitable for laser propagation. These are the 3-5 |lm short wave infrared (SWIR) and the 
8-12 |xm long wave infrared (LWIR) bands. Now that operation at 3 p// with a 1 MW output 
has been demonstrated, operation in the LWIR band at a wavelength about 10 pm, using the 
same undulator, is explored. 

To maintain resonance at A, = 10 pm, the beam energy must be reduced to 54.8 MeV. 
The dimensionless current and optical mode waist are now j = 13 and = .66 mm, 
respectively. If we lay a curve corresponding to a 1 MW output with the new parameters on 
top of Figure 17, we find that I ~ 30 and Q ~ 4.5. Therefore, if the resonator optics can be 
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made to have the same outcoupling at the LWIR wavelengths, the same undulator and 
resonator cavity could be used to operate in both bands. 

The simulation results equivalent to Figures 18, 19, and 20 reveal only a few 
characteristics different than for the A. = 'i\un case. There are four principal longitudinal 
modes excited by the micropulses. Therefore, there is a modulation on the complex field 
envelope. This results in a startup time increase for the FEL of a factor of three. These 
features, while interesting, do not affect the engineering and feasibility results of concern in 
this analysis. 

Despite the apparent wide bandwidth in operating wavelengths, there are many 
repercussions which must be considered. First, to reduce the beam energy by about one-half, 
the accelerator E is reduced in half. Therefore, the energy storage in each accelerator cavity is 
reduced by a factor of four. This means that Ai / lU increases by a factor of four resulting in 
an equal fractional increase in the rf field distortions caused by an accelerating micropulse. 
This will likely be unsatisfactory for proper beam transport in RT cavities. SC cavities will be 
more resistant to this effect. 

A second consideration is the optical mode size with respect to the undulator gap 
distance. In the event of perfect alignment, the undulator will not penetrate into the Gaussian 
profile of the optical energy. However, if the alignment is off by as little as 1 mm. the 
undulator end magnets will penetrate into space normally occupied by several kilowatts of 
power. This can be compensated for by increasing the undulator gap distance. However, as a 
result of the corresponding changes in the undulator parameters, the undulator will be longer, 
have an increased sensitivity to electron beam quality, and the Q value of the optical resonator 
will increase. This will increase the power density incident on the mirrors for the same mirror 
separation. Therefore, if operation in the LWIR is desired, this tolerance must be considered 
when finalizing the undulator and resonator design. 

D. RESONATOR OPTICS 

One of the principal disadvantages of the FEL for use as a laser weapon is that the 
optical waist is very small. Therefore, the power density incident on the resonator optics is 
large. Figine 21 shows the average power density incident on the resonator mirrors as a 
function of mirror separation for the undulator with A.,, = 3 cm and N = 16. With an 18 meter 
mirror separation, the mirror power density is 70 kw/cm ‘ at >\. = 3 pm and 20 kW/cm ' at 
X= 10 pm. Modem optical components can support a power density in excess of 100 
kW/cm- [51]. 
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power density vs resonator length 




Figure 21 - Power density on the resonator mirrors for operation at A. = 3 pm and 
X = to pm ( * - denotes allowable mirror placements consistent with the rf 
frequency choices in Table 7) 



In the shipboard environment, space and volume allocation is of particular concern to the 
ship design engineers. The mirror power density limitations require that the mirrors be about 
18 m apart. However, it may be possible to reduce the resonator length with mirrors that 
exceed the 100 kW/cm^ performance level. While this is a considerable length witliin the ship, 
the actual volume of the resonator is quite small since the optical mode can be contained 
within a cross-section of about 20 cm. The hardware for tlie active stabilization of tlie optics to 
maintain resonator alignment during ship vibrations, flexing, and torsional rotations will 
constitute the bulk of the resonator volume. Therefore, the only significant amount of space 
that will be needed will be at the mirror locations themselves. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



The PEL is a promising candidate for a possible high-average power shipboard laser 
weapon. Unlike other candidates, such as chemical lasers, the PEL is electrically powered. 
An additional electrical power distribution bus would be required aboard a ship if an PEL 
weapon was to be deployed. 

An PEL with a 1 MW laser output would require 6-10 MW of electrical power from the 
prime power distribution system. A gas turbine naval engineering plant could be modified to 
supply this power. Power from a propulsion turbine could be diverted to power the dedicated 
PEL distribution bus. The added distribution system would increase the ship’s displacement by 
22 tons and require 22 m^ of the ship’s volume to implement The power supply that powers 
the rf sources for the linear accelerator were rated at an output voltage of 100 kV and an 
output power of 10 MW. This power supply is only 82% efficient and requires a large 
external cooling capability. The efficiency can be improved, however, with the development of 
alternative rf sources which operate at half of this voltage. The development of these alternate 
rf power supplies is recommended. 

The linear accelerator driving the PEL would employ energy recovery to increase the 
power efficiency and reduce the operating radiation levels. Consequently, the electron beam 
has a large average current. Purthermore, the accelerator would operate at lower rf frequencies 
than most linear accelerators. Therefore, the electron pulse lengths are longer than those 
commonly used in most PELs. The high-average current and long pulse lengths reduce the 
number of undulator periods required to extract the 1 MW output power from the electron 
beam. Therefore, the PEL has a low sensitivity to electron beam quality. 

The PEL examined in this thesis requires further development before it can be built. 
Some of the technologies that require further development are: energy recovery employment, 
control system requirements, and the accelerator electron injector output capacity. Additional 
experiments are needed to determine the limitations of energy recovery including rf phase and 
amplitude stability requirements as a function of PEL undulator energy extraction efficiency. 
Control systems are required to regulate the input power to rf sources and the rf input to 
accelerator components. Additional work is necessary to adequately control the electron beam 
transport instabilities that will result from high-average currents. Pinally, electron injector 
technology can currently produce average currents of tens of milliamps. Advancements are 
needed for the shipboard PEL weapon. Punher development of these FT;L technologies is 
recommended. 
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APPENDIX - POWER SUPPLY CALCULATIONS 



This appendix outlines the procedure used to estimate the weight and volume of the PEL 
power supplies in Chapter IV. Only commercially available, off-the-shelf components were 
considered in the analysis. No attempt was made to determine the savings that may be 
possible with a product development program. Therefore, the results are a conservative 
estimate. The components selected in the analysis are meant to be representative components. 
It should not be assumed that these components are recommended for procurement. The 
voltage and current ratings for all of the components used in the analysis were derated by the 
Navy derating factor of 0.7. The dc link and filter components were not considered. 

A. 100 KV, 3 PHASE, PHASE-CONTROLLED RECTIFIER 

The functional circuit diagram for a 3 phase, phase-controlled rectifier is shown in Figure 

22 . 




Figure 22 - Functional circuit diagram for a 3 phase, phase-controlled rectifier 

Commercially available silicon conu^olled rectifiers (SCRs) are not rated to support a 100 kV 
blocking voltage. Several SCRs must be connected in series to accomplish the function of 
each SCR in Figure 22. Each of these series-connected SCRs must be protected with a 
snubber circuit and balancing resistor. The series connection of the SCRs is illustrated in 
Figure 23. Forced-commutation circuits for the SCRs are not included. It is assumed that the 
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Figure 23 - Series-connected SCRs with protective balancing resistors and snubbers 

blocking voltage is sufficient to rapidly reverse bias the SCRs. 

1. SCR Selection 

The rectifier has an output dc voltage, of 



where is the peak line-to-line input ac voltage and a is the SCR firing angle. Since the 
output dc voltage is reduced by only 13.4% when a is adjusted from 0-30®, V„ = 105 kV is 
chosen so that = 100 kV for a = 0°. To supply 10 MW of power, the output dc current, 
must be 100 A. The peak reverse blocking voltage for the SCR, V„, is V„ and the 

conduction current, /, is 100 A. The derated values of V„ and / are = 150 kV and = 

143 A, respectively. 

Using Ref. [61], the SCR model T8K7 was chosen for the SCR component. The ratings 
of this SCR are listed in Table 14. The value Vrum is the repetitive reverse blocking voltage, 
h{av) is the peak average conduction current, Igt is the gate current, is the leakage 

airrent, is the peak repetitive surge current, and is the turn-on delay time of the SCR. 



^dc = - cos a = 100 k V. 



(A.1) 



n 
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^RRM 


iav) 


hr 


^RRM 


(kV) 


(A) 


(mA) 


(mA) 


4.5 


350 


200 


75 



hsM 


Id 


Weight 


Volume 
(cm -^) 


(A) 


|Asec 


(kg) 


1500 


2 


0.227 


72.9 



Table 14 - Rating values for SCR model T8K7 



Allowing for the failure of one SCR, 35 SCRs are required to construct each of the six 
branches in Figure 22. Therefore, a total of 210 SCRs are needed. The total SCR conduction 
losses are 19 kW. 



2. Balancing Resistors 

Balancing resistors are required to ensure that aU series-connected SCRs block the same 
voltage. Otherwise, variations in leakage currents could cause an excessive voltage to build up 
across one SCR. The actual blocking voltage for each SCR is V'woci = 3.1 kV. The value of 
Rf, is calculated as 

R = = 4 kQ (A.2) 

10 1 ^^^ 



The required power rating of this resistor is 



V 

— = 4.8kW 



(A.3) 



where Vj = 4.5 kV is the derated voltage rating of the SCRs. Multiple resistors will have to 
be placed in series to meet this power rating. From Ref [62], a suitable resistor is the type 
270. The ratings of this resistor are listed in Table 15. The actual power dissipation of this 
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Power 


Weight 


Volume 


(kW) 


(kg) 


(cm 3) 


1 


0.28 


1610 



Table 15 - Rating values for the type 270 resistor 

resistor is Vbiockl^b = 2.4kW . One balancing resistor is needed for each series-connected SCR 
and five type 270 resistors are required per balancing resistor. Therefore, a total of 1050 
resistors are needed with a total power dissipation of 500 kW. 

3. Snubber Components 

In this rectifier, the snubber provides voltage and current protection for series-connected 
SCRs during firing and turn-off. If one or more SCRs latch slower than the rest of the SCRs 
in a branch, an excessive voltage could be applied across these SCRs. The snubber capacitor 
provides the charge to the conducting SCRs during the period when these slower SCRs still 
have not latched. When the SCRs are first conducting, though, the capacitor would be 
effectively shorted causing a large current surge. The snubber resistor is placed in series with 
the capacitor to limit the discharge. 

There is no definitive way to analytically determine rectifier snubber component values. 
However, a conservative methodology can be used to estimate these values with a design 
margin. The following method describes how the snubber component values were determined. 

The latching delay time for the T8K7 SCR is 2 |lsec. Assuming the snubber capacitor 
must supply the full conduction, the minimum capacitance necessary to store the required 
charge is 

Cs^in = = 0.09 \lF . (A.4) 

* block 

The snubber resistor must limit the capacitor current to within the surge current rating of the 
SCR. The derated surge current rating of the SCR, I^sm . is 1050 A. The minimum snubber 
resistance is calculated by 
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= 5Q 



(A.5) 



R 



s jTun 



^dsm ~ ^d 



Assuming that the SCR switches instantaneously, the power rating of the snubber resistor is 
required to be 



^'r, - ^rms ~ 



R. 



K,CJa 



(A.6) 



where fac is the rectifier input ac frequency. Equations A.4, A.5, and A.6 include a margin of 
conservatism. The resistor power losses are independent of the resistance value but are directly 
proportional to the ac bus frequency. To minimize the power losses while maintaining the 
effectiveness of the snubber, the values of R^ = 5 Q. and Cj = 0.2 |iF are required. The 
required voltage rating of the capacitor is 4.5 kV, and the resistor is to have a power rating of 
3.2 kW. The value fac = 800 Hz is assumed. 

From Ref. [63], a suitable capacitor design was selected. The capacitor ratings are listed 
in Table 16. 



Voltage 


Capacitance 


Weight 


Volume 
(cm ■’) 


(kV) 


(HF) 


(kg) 


4.5 


0.2 


0.43 


400 



Table 16 - Rating values for a 0.2 pF, 4.5 kV capacitor 

There are 210 required snubber capacitors. Therefore, a total of 210 capacitors are required. 

The type 270 was selected for the snubber resistor. Four type 270 resistors are needed 
for each snubber. Therefore, a total of 840 resistors are required. The power loss of each 
snubber resistor is Vbiock^sfac ~ l-54kW, so that the total power dissipation by all of the 
snubbers is 1290 kW. 
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4. SCR Gating Power Supply 

A current pulse applied to the gate of an SCR is required in order for it to fire. The 
T8K7 SCR requires a current pulse of 200 mA and dissipates 5 W at the gate in order to fire. 
Since there are 35 SCRs in each of the six rectifier branches, a total of 7 A and 175 W is 
required to fire all of the SCRs in a branch. The derated values of required gate power supply 
current and power are 10 A and 250 W, respectively. 

From Ref. [64], Acopian power supply model A28H1400 was selected. The ratings of 
this power supply are listed in Table 17. 



Voltage 


Current 


Power 


Weight 


Volume 


(kV) 


(A) 


(W) 


(kg) 


(cm^) 


28 


12 


336 


11.8 


9900 



Table 17 - Rating values for Acopian power supply model A28H1400 

A total of six power supplies would be required. 

B. VOLTAGE MULTIPLICATION CIRCUIT 

The fundamental circuit diagram for a 500 kV voltage multiplier is shown in Figure 24. 
Figure 24 is an 8x voltage multiplia. An increased multiplication factor, F^, is obtained by 
cascading additional stages. To obtain a 500 kV output with an input line voltage of 5 kV 
(rms), a multiplication factor of F^ = 72 is required. It is important to know that if any of 
these stages fail, the voltage multiplication factor is reduced by two. The next stage in the 
cascade must be rated to assume the load of the failed stage or it may also fail. Therefore, 
each diode and capacitor must be designed to the limiting position in the circuit. The diodes 
must be designed for a reverse blocking voltage of 2£^ and a current of The capacitors 

must be rated for a voltage of 3Ep . The derated peak surge current, to which aU of the circuit 
components must be rated, is given by 

idsn, =^facCEp=22S.5A. (A.9) 
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Figure 24 - Functional diagram for a voltage multiplier circuit 



1. Capacitor selection 

The effective output capacitance. Cefj , of the voltage multiplier needed to supply a 
current Ijc is given by 



C 



/, 



dc 



<’// 



foe 



(A.7) 



where is the allowed output voltage ripple. The value of the circuit is 



Cejj 



= 4 C 



^-1 

F^F,-2) ' 



(A.8) 



where C is the capacitance of each capacitor. Given values of = 800 Hz. = 500 mA. 
Fy = 72, and arbiU'arily choosing AV''„ = 5 kV. then each capacitor in Figure 24 must have a 
capacitance C = 2.25 pF. The derated voltage rating of the capacitors is 30.3 kV. The total 
number of these capacitors is 1.5F,, - 2 = 106 capacitors. 

Using Ref. [63]. capacitor model 30F1377 was selected for the capacitor component. 
The ratings of this capacitor are listed in Table 18. 
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voltage 


Capacitance 


Weight 


Volume 


(kV) 


(ItF) 


(kg) 


(cm^) 


30 


2.37 


19 


11500 



Table 18 - Rating values for capacitor model 30F1377 

A total of 106 capacitors is required. 

2. Diode selection 

Each diode in Figure 24 must be rated for the limiting current and voltage values in the 
circuit. The required diode ratings are = F^Ij^/0.1 = 51.4 A, = 20.2 kV, and - 
228.5 A. Using Ref. [61], diode R4051070 was selected. The ratings of this diode are listed 
in Table 19. 



^RRM 


(av) 


hsM 


hm 


Weight 


Volume 


(kV) 


(A) 


(A) 


(mA) 


(kg) 


(cm^) 


1 


70 


1200 


0.1 


0.018 


3.18 



Table 19 - Rating values for diode R4051070 

It is necessary to place several of these diodes in series to accomplish the function of each 
diode in Figure 24. Allowing for a single failure within the diode branch, 22 diodes are 
required per branch. There are F^, diode branches so that a total of 1584 R4051070 diodes are 
required. The total conduction losses are 50 kW. 
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3. Balancing resistor selection 

The peak blocking voltage of each diode is Vhiock = 2£^/22 = 643 V. Using equation 
A.2, the required value for the balancing resistor is Ri, = 643 kfi. The required power rating 
for this resistor is VjiRf, = 1.3 W. Using Ref. [65], resistor 620KW-2-ND was selected. The 
ratings of this resistor are listed in Table 20. 



Power 


Weight 


Volume 


(W) 


(kg) 


(cm^) 


2 


0.01 


0.4 



Table 20 - Rating values for resistor 620KW-2-ND 

One of these resistors is required for each of the R4051070 diodes. Each of these resistors 
dissipates Vhiocklf^b =0.64W. Therefore, 1584 resistors are required with a total power 
dissipation of 1 kW. 

C. SUMMARY 

The total weight, volume, and power losses for the 100 kV rectifier and 500 kV voltage 
multiplier power supplies are listed in Table 21. 





100 kV 


500 kV 


units 


Weight 


1500 


4100 


kg 


Volume 


8 


3 


m'^ 


power losses 


1800 


50 


kW 



Table 21 - Power supply summary 
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The weight estimate was made by multiplying the component totals by a factor of two. 
Power supplies that operate above about 40 kV are required to have enclosures which are filled 
with high voltage oil, SF g, or some other means of protection between components. Due to the 
large amount of heat losses, the 10 MW power supply is assumed to be filled with oil which 
can also be used as a cooling medium. The 500 kV power supply will need high voltage 
shielding to isolate it from the environment. The volume estimates include a stacking factor of 
2.5. This factor accounts for the high voltage shielding and unused space within the power 
supply after the components are assembled. 

The power efficiency of the 10 MW power supply is about 82%. While this is 
acceptable for desktop size electronics, the 1.8 MW of power losses place an added design 
burden of forced cooling for the power supply. The efficiency of the power supply can be 
increased significantly with the development of alternative accelerator rf sources, as described 
in Chapter III, which operate at substantially lower voltages. 
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